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ABSTRACT 


Relativistic wave functions are derived for the electron in a potential 
field which is constant within the nucleus and Coulombian outside. The 
relative values of the cut-off wave-functions and the pure Coulomb 
wave-functions are calculated both for the discrete and continuous 
spectrum. The changes thereby produced in the probability of orbital 
electron-capture, in the continuous Beta-spectrum and in the ratio of 
probabilities of K-capture to positron emission are examined. It is 
concluded that for the heavier nuclei the probability of K-capture and the 
ratio of the probabilities of K-capture to positron emission are appreciably 
affected. The modifications produced in the continuous spectrum are 
greatest for spin change equal to n where n(n>1) is the degree of 
forbiddenness of the transition. In this case the coefficients of the 
squares and products of matrix elements given in the correction factors 
of Konopinski and Uhlenbeck are appreciably changed. Therefore for 
given values of the matrix elements the spectrum is modified. 


§1. INTRODUCTION 


In a previous paper by Malcolm and Strachan (1951), an investigation was 
made into changes produced by a cut-off field in the relativistic wave 
functions of the bound electron. Results were computed only for the 
case of the K-electron. The calculations have since been extended by the 
present author and, in table 1, for the K, Ly and Ly, electrons, and for 
different values of the atomic number Z, are given the values of the ratios 
LPR, gly? and y=C(P+Ig PMP P+I 9°} evaluated at the 
nuclear radius where f, g, f°’, g° are the radial wave functions for the 
cut-off field and the pure Coulomb field respectively. 

From examination of table I we find, as we would expect, that | gz, |? 
is affected slightly less than |g,|?. However, it is interesting to note 
that, although |f, |? and |f;,|? are very much affected by the cut-off 
field, | f,,,[?, contrary to expectations, is not, its decrease being just 
slightly less than that of |gz,|?. We can therefore conclude that the 
cut-off field has most effect on K-capture and that its influence on orbital 
decay becomes progressively less for succeeding orbits. For L-capture no 
comparison with experimental results is made here. However, it is to be 
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noted that, for the heavier nuclei, the effect of the cut-off field on the 
probability of L-capture will be greater than the order of magnitude of the 
change we would expect from screening (Slater 1930). The value of the 
ratio of L to K-capture for any nucleus will be changed only slightly by 
use of the cut-off functions for, as one can see from table 1, their 
probabilities have approximately the same decrease. 

The components of the Dirac wave function which are affected most 
do not occur in the calculation of orbital electron capture. Those 
components are the third and fourth when « is negative and the first and 


Table 1 


Z, A 40,89 | 48,107 | 50,113 | 60,146 | 70,170 | 80, 204 


| faci be 0-400 0-383 | 0-376 | 0-326 
[fc,/f%, 0-418 0-405 | 0-387 | 0-364 | 0-333 
Lf taal Sn 0-980 0-968 | 0-952 | 0-931 | 0-902 
lglg lk : 0-930 | 0-904 | 0-863 
| 91,/9%, be 0-976 0-963 | 0-944 | 0-920 | 0-887 
{ 9a/9n | 0-451 0-433 | 0-414 | 0-390 | 0-358 
Xx 0-903 | 0-868 | 0-812 
Xt, 0-964 0-943 | 0-916 | 0-875 | 0-834 
Xn 0-966 0-947 | 0-922 | 0-884 | 0-846 


second when « is positive ; «——(l+-1) or/ according as the total angular 
momentum quantum number j is equal to /+4, /—} respectively. 
However, for the continuum all components of the electron wave-function 
may play an important part in the determination of the energy distribution 
of forbidden transitions. We therefore determine the cut-off wave 
functions of the electron for the continuous spectrum. 


§2. WaAvE-FuNcTIONS 


Energy, momentum and length are measured in units mc?, me, h/me 
respectively. The equations to be solved are 


d 
4 + (11) fir (WY) p=. Re era (6, 
d. 
t+) gt—WW4+1—T)f=0,, . . . . @) 
where View ae wae (3) 
—Zalr (R<r: Region II), 


a=e/he, 


Here the signs are in accordance with Rose’s (1937) definition of the 
spherical harmonics which give the angular dependence of the wave 
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functions. In Region I the solutions which behave suitably at r=0 are 
formally the same as those obtained for the discrete spectrum, viz., 


f=A,fr=Ay™ J iou(p?), 


"A (4) 
9= Ay =F sa (Ps 
La. 
where p=1/(W"?—1) and WW =. Pe ork ao} 
The general solutions in Region II can be written as 

f=A(futkFn), | (6) 

g=A (9 +hGyy), 
where fq, 91, Fy, Gy, are defined as follows : 
fr | =F W)P(Ripry 
In {Lexp(—ipr+%) (y—8)F(1+y—8, 1+ 2y, 2tpr)Fe.c.}, 5e i) 
Fy | =F W)*(2ipr)-?y—8) 
Gry {L[exp(—ipr-+in)|F(1—y—d, 1—2y, 2ipr)-e.c.} eet ama. § (3) 


where c.c. stands for the complex conjugate and F is used for the hyper- 
geometric function ,/’, in the notation of Pochammer and Barnes (Watson 


1944). Also 
: 7 
pav(Wi1), b= AO, a payee). (0) 
and exp (2i7)=(—«+8’)/(y—8). 


We note that, apart from the normalization constant, f,; and g;, are the 
solutions as given by Rose of the Dirac wave-equation for a pure Coulomb 
field. When normalized our wave functions in Region II may be written 


; f=(NS)(fut+kFu), 
g=(NoS)(Gu +k), 


- where (NS) is the normalization for the cut-off potential and Ng is the 
usual Coulomb normalization factor. Since (No¢fjr), (Negi) are the 
normalized pure Coulomb wave-functions, say f*, g° respectively, we obtain 


falps (14452), 


(10) 


Gg (11) 
=(9°)S (1 bot) : 
g=(9°) ( Tahoe 
The boundary condition at r=, the nuclear radius, gives 
need Irfu—hGn ) (12 
S Bae Gu) Rr ) 
322 


e 
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Using this substitution for k in (11) gives for the value of our wave- 
function at the nuclear radius 


(13) 


where Q=(fi 9ul9 fe T=(fi Gulgr Fu)r- at fh ars ke 


Hereafter the suffices -+ affixed to Q and 7’ will mean the values of these 
quantities for « positive and x negative, respectively. Analytical 
expressions for Q,, and 7’, are given in Appendix L 


§ 3. NORMALIZATION 


For normalization of the continuous wave-function we arrange that 
[8 dG w dew thawte wI=N3W'—W), . . « (15) 
0 . 
where 8(z) is the Dirac delta function. Along with eqns. (1) and (2) 
this requires 
Whee va : : 
N= | dW’ lim Wow 95 w lew Se, Ix, w) aes (16) 
w—Aw r—>o 
where A of equation (7) satisfies A=1/,/N’. 
Evaluation of the integral (16) gives 


(N’Jeuton=(I+ED\(L+-K"D"\(N")coutomy + + + + (17) 
P+y—s) (1-2 
where De eae ay < joa! We hee 
Therefore, for the cut-off field, A is given by 
A=N(S=N f1+-kD)ULED 3) ee 
where N¢ is the pure Coulomb normalization factor. Expanding by the 
binomial theorem and dropping the term | kD |?/2 which will be very small 


we have S=[1—AkD))\ 2 Oe ee 


where Z{ kD} stands for the real part of the argument. 

The above method gives the normalization factor for solutions normalized 
per unit energy interval. To obtain solutions normalized to one particle 
in a sphere of unit radius one must multiply the above solutions by (zp/W)1?. 
We shall assume that this has been done and that the extra factor has 
been absorbed into the pure Coulomb normalization factor Ng. Formulae 
fork, and &{kD} are given in the Appendix. Calculations were carried 
out for |«|>1, but the resulting change in the normalization factor was 
found to be negligible. Computation of the values of [1—A{kD}}?? (over 
the appropriate energy ranges for the disintegration) was carried out for 
Z=44, 51, 70 and 84 and | « |=1. It was found that even for | «|=1 and 


high Z the change in the normalization factor squared is very small, being 
at most about 1%. 
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§4. Resutts or Computation 


Using Rose’s suffix notation for f and g, we obtain from eqns. (13) 
the following relations between the cut-off’ wave functions and the 
Coulomb wave-functions evaluated at r=R :— 

For « negative, 


(fila==CIff)n where Of—S, (1 - =) 


1—T_ 
(21) 
LSA g ae 
(W)r=CU9Dr Where Cl=S, i= 1_7T-1]° 
For « positive, 
' [= 
area O42. fo 725) x where CH Spier (1- l Tr) ) 
eee 
7 eee (22) 
(9-1-1) r=C” 1_-1(97:1_1) rz Where Ce 82-4 l= = : 
IT 
Also S=(1-Z2(k_D}), S_)-1=(1—A{keD}) . (28) 


where k_, k, and D are defined by eqns. (A1), (A2) and (18) respectively. 

For four values of Z, calculations were made of the values of 
(Cf)?, (C4£,_,)?, (C4)? and (C”_,_,)? for 1=0, 1, 2. The values of Z chosen 
give a reasonable distribution in Z; other reasons for this choice will be 
discussed below. In table 2 the results of these calculations are given in 
full for the case Z=84. For the other values of Z we quote in table 3 
only the mean values of the above quantities over the energy range 
considered, since in all cases it was found that the C/ and C” are almost 
independent of the energy. 


Table 2 


is (OL)? | (O2,)2 | (Of)? | (O%,)2 | (C22 | (CL,)2| (CN? | (CL,) 


0:3274 | 0-3529 | 0-6170 | 0-6350 | 0-8796 | 0-8969 | 0-9882 | 0-9892 
0-3295 | 0:3534 | 0-6189 | 0-6367 | 0-8776 | 0-8934 | 0-9880 | 0-9890 
0-3301 | 0:3536 | 0-6210 | 0-6385 | 0-8751 | 0-8911 | 0-9878 | 0-9889 
0-3307 | 0-3538 | 0-6228 | 0-6401 | 0-8727 | 0-8883 | 0-9877 | 0-9887 
0-3313 | 0-3542 | 0-6247 | 0-6419 | 0-8704 | 0-8856 | 0-9875 | 0:9886 
0-3319 | 0-3545 | 0-6265 | 0-6437 | 0-8682 | 0-8830 | 0-9874 | 0-9884 
0:3326 | 0-3549 | 0-6283 | 0:6454 | 0-8662 | 0-8805 | 0-9872 | 0-9883 
0-3333 | 0-35538 | 0-6302 | 0-6471 | 0-8642 | 0-8782 | 0-9870 | 0-9881 


From examination of tables 2 and 3 we can therefore draw the following 
conclusions : 

(a) f,, the small component when « is negative, and g_; 1, the large 
component when « is positive, are very much affected by the deviation 
from the Coulomb field. This effect is appreciable even for a nucleus 
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whose value of Z is as small as 43 (see table 3). The other components 
are much less affected, but it is interesting to notice that for Z as high as 
83 (see table 2) even | Jo |? has been reduced by about 12%. 


Table 3 


Decaying| (og | (02,2 | (Of? | (C25)? | (OR? | (CZ) | (CD® | (CZ)? 


nucleus 


al Bs 0-4200 | 0-44.94 | 0-6370 | 0-6620 | 0-9763 | 0-9799 | 0-9969 | 0-9973 


1258 | 0-4314 | 0:4575 | 0-6530 | 0-6718 | 0-9641 | 0-9952 | 0-9681 | 0:9957 


170Tm | 0-3761 | 0:4059 | 0-6320 | 0-6512 | 0-9261 | 0-9359 | 0-9916 | 0-9925 


nn rE 


(b) As one would expect, the effect of the cut-off field on the components 
of the wave-functions decreases as the angular momentum increases. 

(c) The values of the C/ and C% as defined by eqns. (21) and (22) are 
almost independent of the energy over the ranges considered. 


§5. MopIFICATION OF THE FORMULAE OF KONOPINSKI AND UHLENBECK 
(1941) 
We first consider the effect which the above results will have on the 
L, M, and N, as defined by Konopinski and Uhlenbeck. Their formulae, 


expressed in terms of the Coulomb functions, become for the cut-off 
function, 


L=u(Z, W)\(C4)? {(9%)?+03(f2 1-1)" be 
Mi=0(Z, W)(Cf) {(F§)? +00(9%-1-1)? he: ee re eS 
N,=0(Z, W) CLO4{figi—asf21197-1-1 he | 
where He eee 
uo(Z, W)= (Fre, w)) ae a= (C1, ,/C))*, a—(Co op 


tg=(CL,_,0%,_,/CfCP). 


Examination of tables 2 and 3 shows that «,, x) and a3 are approximately 
equal to unity and that they vary only very little with energy. Therefore, 
since the C’’ and C/ have also only a slight energy dependence, the major 
effect of the cut-off field on the L, M, and N is to multiply them by a 
factor which is approximately constant over the energy range. : 

We represent the forbidden correction factors given explicitly by 
Konopinski and Uhlenbeck in the following manner : 


C, = ZA;(A; shape) ".7°.09. 0 eee ao) 
where the index n denotes the order of forbiddenness and the Y represents 


. . y 
_ a summation over the various squares and products (A,”) of the matrix 
elements which occur. In (25) we have not used an index to denote the 
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type of interaction, since, whenever we need to particularize in the ae 
this will be made clear. 
For allowed transitions YA, (A,° shape) reduces to one term, say 


r 
A,° (A,° shape), where (A,° shape)=L,. For first forbidden transitions 
the (A," shape) can contain only sums of Lo, L,, Ny and M,. Also for 
heavy nuclei the M are the largest terms ; for example, 


Za\2 Za 
Me. bre NG (x) yar Il, 


If we let ¢ denote the neutrino energy, we can write the different types of 
(Af shape) that occur as follows : 

(a) (A,1 shape)=M,+2qN,+terms in Lo, L,, 

(b) (A,1 shape)=—N,+3qLZ,, 

(c) *(A,1 shape)=L,, 

(d) (A,1 shape)=5¢?L,+3L,=« factor. 

We can therefore deduce that (A,° shape) and (A,1 shape) type (c) 
will not be affected appreciably by our modification since they depend ona 


mie 


Wunit=mc?) 


1-0 2-0 3-0 
RaE. Curve (a) is the a factor curve unmodified ; 
curve (0) is the modified « factor curve. 


single term which, when corrected for the influence of the cut-off field, is 
virtually multiplied by a factor independent of the energy. However, 
A, shape type (d) contains two terms one of which (¢°L) is dominant at 
low energies and the other (?1,) dominant at high energies. Therefore, 
if the reductions in L, and L, are appreciably different, (A,1 shape) for this 
type of transition will be modified. This is true for Z=84 and for Z=71 
although for the latter the effect will be much smaller. Since (C{)?> (C$)? 
the corrected (A,! shape) will have a relatively greater depression at the 
low energy end as compared to the high energy region. This is shown in 
fig. 1 where for Z=84 we have drawn the modified and unmodified « 
factor curves. The above argument applies in a lesser degree to (A,1 
shape) types (a) and (b). Here, although the leading terms are largest at 
- all energies, the effect of the second terms may be appreciable at low 
energies or more especially at high values of q. We would therefore 
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expect that, for moderate values of the end point energy, the effect of the 
cut-off field on the actual shapes of the (A,! shape) types (a) and (0) will be 
intermediate between that produced on types (c) and (d); however, 
the (A,1 shape) types (a) and (b) will be reduced by approximately constant 
multipliers which deviate more from unity than those for types (c) and (d). 

The modification of the (A, shape) type (d) depends on the changes 
produced in the ZL which are the least affected by the cut-off field. 
Therefore for the maximum effect of the cut-off field we must examine 
the second forbidden transitions or higher. For second forbidden 
transitions the (A,2 shape) can contain only sums of Lo, Ly, Ly, No, Ni, 


Fig. 2 


100 (2) 


50 
(2) 


W(unit = mc?) 
1-0 2:0 3-0 
RaE. Curve (a) is (2|7'j|? shape) unmodified ; 
. 


J 
curve (b) is the modified curve. 


M, and M,. We write them as follows, only the largest terms being 
given explicitly when much smaller terms also occur. 
(a) (A? shape)=4(6q'L)+2¢?L,+15L,)=D, factor, 
(6) (A? shape)=5(@M,)+9M,+....), 
(c) (A,? shape)=« factor, 
(4) (A,2 shape)=}(—4¢?N)—3N +... -), 

(e) (A? shape)=(M,+2qN + ....). 
We have already considered the change in types (c) and (e). Type (a) 
is somewhat similar to the « factor and for Z as low as 44 and 51 the 
modification produced in it is negligible. The maximum influence of 
the cut-off field on the (A,? shape) is exhibited by terms of types (6) 
and (d), especially type (6). These are the coefficients which occur 
for n=2 and spin change of 2 (i.e. unfavourable parity change in the 
case of Gamow-—Teller interactions). In types (6) and (d) the (A,2 shape) 
are dominated by q?M,, —4q?No respectively at low energies and by 
9M,, —3N, respectively at high energies. Since C,f and C\f are 
appreciably different and O,/>C,/ it is evident that these (A? shape) 
_ will be affected in the same manner as the « factor though to a much 
greater extent, i.e. there will again be a relatively greater depression at 
the low energy end as compared with the high energy end and the change 
will be the more obvious the higher the maximum energy W,. This is 
illustrated in figs. 2, 3, 4 and 5, where, for Z=51 and 84, the (A,? shape) 


Effect of a Finite Nucleus on Beta-Decay 1019 
Fig. 3 
(2) 


5-0 


Wtunit = 707) 
1-0 2-0 5:0 
Rak. Curve (a) is (2(7,Ajj;+ ¢.c.) shape) unmodified ; 


curve (b) is the modified curve. 


Fig. 4 


W (unit= 7c?) 


1-0 2:0 3:0 40 5:0 
225Sn, Curve (a) is (2|7'j|? shape) unmodified ; 


curve (5) is the modified curve. 


Fig. 5 


W (unit =2c) 
5:0 


1-0 2:0 3:0 4:0 5:0 
1258n, Curve (a) is (X(T yAjj+c.c.) shape) unmodified ; 


ij 
curve (6) is the modified curve. 
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types (b) and (d) have been drawn for the modified and unmodified 
functions. The type of change of shape shown there would seem to hold 
good for any (A,” shape) which is appreciably modified by the cut-off 
field. 

However, we can draw no general conclusion as to the manner of the 
change in shape of the total correction factors which include the 
(A,2 shape) types (b) and (d) since the coefficient of the (A,* shape) 
type (d) may be either positive or negative. The resultant change of 
shape will therefore depend on the size of the matrix elements. However, 
we can conclude with certainty that, for a fixed set of matrix elements, — 
the second forbidden correction factor with spin change 2 should give 
very appreciably different curves for the modified and unmodified 
functions, as well as a change in the ft value. This statement is also true 
for any forbidden transition of degree n (n>1) and spin change n. 
It should be noted that even in the cases where the shapes are not 
modified the ft values will be changed. 

While the work described above was being used for comparison with 
experimental spectra (see § 6), a letter by Rose and Holmes (1951) was 
published which gave the results of their consideration of the same 
problem. The conclusions they reach are in the main in agreement with 
those stated above. However, we can see no reason for the general 
statement made at the top of p. 91, column 2, that ‘a relatively greater 
depression of the spectrum is to be expected at the low energy end as 
compared with the high energy region’. For the orders of forbiddenness 
which Rose and Holmes consider this seems to the author, as stated 
above, to depend on whether the cross term is negative or positive and 
on the order of magnitude of the different matrix elements. 

Rose and Holmes used approximately a uniform charge distribution 
within the nucleus (actually they include a term representing charge 
distortion but find their results insensitive to this) in place of the charged 
shell at the nuclear radius considered here. This means that within the 
nuclear radius their potential well is approximately half as deep again 
as that used here. This difference affects some results considerably. 
For example, for Z=83 (RaE) Rose and Holmes calculate that 
(Co’)?==0-47 while for the constant potential inside the nucleus we find 
(Co/)’==0-33 ; similarly for (C/)? Rose and Holmes calculate (C,/)?==0-8 
against our value of (C,/)’=0-62. (The approximation signs have been 
used because a rough mean value over the energy range has been quoted.) 
However, since the ratio (C/)"/(O,/) is approximately the same in both 
cases, little change in the shape of the spectrum for the two cases will 
result though the ft values will be different. Rose and Holmes point out 
that their modification may only be important for Z>~60. However, 
because of the shallower well employed above, lighter nuclei will be 
affected more in our calculations, and this is one reason those calculations 
were made for Z with as low a value as 44. 


ee x. Sh Ue CU 
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§6. APPLICATION TO SPECTRA 


We have already shown that the cut-off field has an important influence 
on certain types of Beta-transitions. The difficulty is to illustrate it 
by application to particular cases. The surest way would be to calculate 
the matrix elements for a particular transition on the basis of some 
nuclear model, although this would be very difficult for some of the 
heavy nuclei with which we are dealing. Another method, which was 
the one decided upon, is to attempt to fit the modified and unmodified 
correction factor curves to correction factor curves derived from experi- 
mental results. This method seemed to be the better because (a) one 
could compare the two best fits so obtained and (b) even although there 
may be little choice between these best fits (as might easily happen since 
we are allowing the matrix elements to vary), the values of the matrix 
elements so derived would illustrate the difference produced by considera- 
tion of the finite nucleus. They could also be compared with theoretically 
derived values of the matrix elements if so desired. 

The values of Z used for the calculations of § 4 were chosen because 
those values of Z belong to the product nuclei of disintegrations with 
reasonably reliable experimental spectra whose shapes have so far, in 
our view, not been adequately explained and which might be explicable 
by applying the modified correction factor for second forbidden transitions 
with unfavourable parity change under Gamow-Teller rules. It was 
decided to use the pure tensor form of the interaction in the fitting because 
evidence indicates that the tensor interaction must be at least part, if not 
all, of the interaction necessary for the explanation of Beta-decay. The 
attempt was therefore made to fit our theoretical correction factors for 
99Tc, 25Sn, 1°°Tm and 2!°RaE to those derived from experiment (see the 
Appendix) ; the correction factors were calculated on the assumption that 
all decays were second forbidden with spin change 2. (Evidence for the 
order of forbiddenness is given below.) 

Using the A,? as unknown constants, the method of least squares was 
used for fitting + A,? (A,? shape) to the experimental curve. Here a 

ue 


difficulty arose which the author has never seen mentioned elsewhere. 
It was found in all cases that, on reduction, the matrix of the equations 
proved to be singular within the accuracy of the calculations ; that is, in 
each case there was virtually a linear relationship between the different 
(A? shape). This is rather surprising, since we are using only three or 
four curves in the fitting. However, it is due to the different (A,” shape) 
having somewhat similar shapes. This singularity of the matrix means, of 
course, that fitting the experimental curve will not determine the ratio of 
the matrix elements uniquely. Whether the shapes of the correction 
factors for the cut-off and pure Coulomb fields are appreciably different or 
not, the ratio of the matrix elements derived from the fitting will certainly 
be much affected. However, comparison of these ratios with each other 
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and with an independent estimate of the matrix elements (say from a single 
particle model) now becomes impossible because of their lack of uniqueness. 

In the forbidden correction factors certain inequalities hold between 
some (A,” shape) ; in particular for the second forbidden factor 


{Z| Ty PHZ|AgP}PETAgtee}. . . . (26) 
ij a aj 


Some authors quote the values of their matrix elements, derived from curve 
fitting, with the subsidiary condition 


4{2|T 4, P}{2| 45 PR={27 yAgtec.? 20 ae oe 
ij yj ij 


which may or may not be true. This condition would, of course, overcome 
the lack of uniqueness in the matrix elements. As an example, in the 
case of RaE, we have also made our fitting subject to (27) rather than (26) 
and quote the ratios of the squares of the matrix elements so derived. 
However, one should not forget that, for fixed values of the matrix elements, 
the actual correction factor curves for the modified and unmodified field 
will be different. 


Technetium ({3Tc) 

The energy spectrum of *°Tc¢ was measured by Taimuty (1951) using 
fairly thin sources which varied in thickness from 0-17 mg/cm? to 
1-3 mg/cm?. He determined the end point energy to be 0-292+-0-003 Mev 
and using this value of W, and a half-life of 2:12 x 10° years he obtained 
an ft value of 2-310" for this decay. Empirically this value would 
make the transition about second forbidden, which is very satisfactory 
from the point of view of applying our modifications. 

Taimuty attempted to fit his spectra by applying second and first 
forbidden correction factors which contained only single matrix elements 
by Gamow-Teller rules. The fits which he obtains are not very satis- 
factory. On the other hand, examination of fig. 6 shows that a very 
good fit has been obtained both for the modified and unmodified Ogp 
correction factor with spin change 2. Hereafter the modified factor will 
be denoted by (Cy7)'.. Of the two, (C7) gives the better fit, the ordinary 
correction factor falling away slightly at low energies. This slight rising 
of the (Czp)’ curve with respect to the Cy, curve at low energies occurs 
also for the other nuclei (except in the case of RaE, where both curves 
are the same). It is interesting to note that this deviation is in the 


same direction as that which would be produced by source thickness 
distortion. 


Tin (*22Sn) 

This nucleus has an ft value of 3-7 108 which suggests that it decays 
by a second forbidden transition. Its energy spectrum has been measured 
by Hayward (1950). From the shape of his Fermi plot Hayward considers 
that the transition is complex and attempts to fit the high energy 
component by the « correction factor. Hayward concludes that the 
fit is satisfactory, but an examination of his figure 2 shows that, in fact, 
the corrected points lie on a smooth curve other than a straight line. 
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For this nucleus also, we have obtained a better fit, with both (Con) and 
Cz, than has hitherto been achieved. Here, however, both Cop and 
(C2r)’ lay within the experimental limits. 


Fig. 6 
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Tc. Dashed curves are the correction factor curves derived from experiment. 


Curve (a) is the best fit for the unmodified C, curve ; curve (b) is the best 
fit for the modified C, curve. 


Fig. 7 
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RaE. Dashed curves are the correction factor curves derived from experiment. 
Curve (a) is the best fit for the unmodified C7 curve ; curve (0) is the best 
fit for the modified C7 curve. 


Radium E 

The first attempt to fit the RaE spectrum was made by Konopinski 
and Uhlenbeck (1941). By use of the C27 factor for spin change 2— 0 
they obtained quite a good fit with the experimental spectrum as 
measured by Flammersfeld, but not such a good fit with the spectrum as 
measured by Neary. Later work, however, by Langer and Price (1949) 
has confirmed the experimental work of Neary. This, therefore, seemed 
an interesting case for the application of our modifications. Figure 7 
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shows the fit obtained with the factors Cy, and (C27)’ (for spin change 
2+ 0) using the results of Langer and Price. Both fits are practically 
the same and could be classed as quite good, though perhaps not really 
as satisfactory as we should like. We obtain for the ratios of 
| A, [/2|7; |? the values 2-9 and 19-7 for the cut-off and pure Coulomb 
ij aj 


U 

fields respectively. Konopinski and Uhlenbeck obtained a value of 
33-6 from their fit. Rose applied his nuclear field modifications (see § 5) 
in an attempt to explain the RaE spectrum. He, also, failed to obtain 
a satisfactory fit for any of the pure invariant forms of the interaction. 
However, he states in his letter that 8. Moszkowski communicated to 
him privately more recent evidence which points to a complex spectrum 
in this case. 


§7. Ratio or K-CapTuRE TO PosITRON EMIssION 


We find that for Z=50, Ly is reduced by a factor 0-94 (a value which 
is almost independent of energy release) for negatron emission. We may 
therefore expect that for positron emission L, will be increased by a . 
factor of approximately 1-06. It was pointed out by Malcolm and 
Strachan (1951) that the experimental value of the ratio of the probabilities 
of K-capture to positron emission for 1°’Cd is 320-+-30, whereas Bouchez, 
de Groot, Nataf and Tolhoek (1950) calculate it to be 360. Taking into 
account only the modifications of the discrete spectrum, our wave 
functions give 343, a value within the experimental limits of accuracy. 
Modification of the wave-function of the positron can only place this 
value more securely within the experimental limits and, if we assume 
that Ly is increased by a factor 1-06, the value 324 replaces 343. As was 
pointed out before, the nucleus 1°’Cd is the only one considered by Bouchez 
et al. which has Z large enough to be affected by our modifications, and 
it is also interesting to note that it is the only nucleus for which they found 
the theoretical ratio of the probabilities of K-capture to positron 
emission greater than the experimental ratio. 


§8. FuRTHER WORK IN PROGRESS 


Very recently, evidence has been accumulating (Horie and Umezawa 
1951, Tolhoek and de Groot 1951, Pursey 1951) in favour of using a 
combination of interactions, and in this connection cut-off wave-functions 


for the electrons are now being applied by the author, especially to the 
work of Smith (1952), 
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BP POND) EX 
Analytical Expressions for Q.., TP ,, ks and kD. 
3 “ad BF W+1\12 
If we let b.(W)= oe A Weal A 
we can derive the following Pies for Q@. and 7’. by expansion of 
the confluent hypergeometric functions occurring in the wave functions, 


We obtain 
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. Le 
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Since y=1/(x?—Z2?), sin 8 is usually small. We have therefore included 

a further term in the expansions of 6 and b’ for use when a and a’ are 

multiplied by sin f (i.e. in the numerators of Q + and his, respectively). 
From eqns. (12) and (14) &., may be written in either of the equivalent 


ENG) OF) lw 


no (FE) (9), (CEE) «= a 
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where (Fu)r=(Fin)a(y—8)(2epR)— “if | 
(Gy) r=(Gia(y—8)(2tpR)- v6: 


To the same order of approximation as was used for Q, we obtain for 
(g11/GS1) zn» (ful Ft) x, zn Tespectively, the following formulae :— 


On : 4 asin B+6b cos B \ 

= UE OS 7 er rae een OD 
= a’ cos a+b’ sin 

On a cos «—6 sin « 
(7) a’ sin B—b’ cos Bf ’ 

hey o- a cos B—6 sin B- 
i a a’ sin «—b’ cos « 
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Except for I'(y—s)/I[(1—y—s5) we have all the explicit formulae 
required for the complete evaluation of the Aik, D} and therefore of 
the normalization factor S. 


where 7.2=y?-+ 


For | «|=1 
ys 3) ae een onl 
Ta a0) a initia aes = y—S) 
iy payee eer, eC 
=(1—y vo (aie 


by application of Taylor’s expansion to the I’ functions. (—8&) is the 
logarithmic derivative of the I’ function and may be expressed as 


z 
MOS Ot Sante)? 


where C=0-57722 (Jahnke and Emde 1945). Therefore 


IB Gp=8 
Popa = 8)(1201—y) (9). 


For |«|>1 one obtains slightly more complicated expressions because 
the I” functions have to be reduced till they are amenable to expansion 
by Taylor’s Theorem (i.e. until the I functions can be expanded in 
powers of |«|—y). Calculations were carried out for |x| >1 but the 


resultant change in the Coulomb normalization factor was found to be 
negligible. 


Effect of a Finite Nucleus on Beta-Decay 1027 


REFERENCES 
Bovucuez, R., pe Groot, 8. R., Natar, R., and Totnork, H. A., 1950, J. Phys. 
Radium, 11, 105. 
Haywarp, R. W., 1950, Phys. Rev. [2], 79, 409. 
Horie, H., and Umrzawa, M., 1951, Phys. Rev. [2], 88, 1253. 
Konopryskt, E. J., and UHLENBECK, G. E., 1941, Phys. Rev. [2], 60, 308. 
Lancer, L. M., and Prics, H. C., 1949, Phys. Rev. [2], 76, 641. 
Matco.rm, I., and Stracuan, C., 1951, Proc. Camb. Phil. Soc., 47, 610. 
Pursty, D. L., 1951, Phil. Mag. [7], 42, 1193. 
Rosz, M. E., 1937, Phys. Rev. [2], 51, 484. 
Rosz, M. E., and Hotmss, D. K., 1951, Phys. Rev. [2], 83, 190. 
SLATER, J. C., 1930, Phys. Rev. [2], 36, 57. 
Situ, A. M., 1952, Thesis, Aberdeen. 
Tarmourty, S. I., 1951, Phys. Rev. [2], 81, 461. 
ToLHokrK, H. A., and pE Groot, 8. R., 1951, Phys. Rev. [2], 84, 150. 
Watson, G.N., 1944, Theory of Bessel Functions (Cambridge : University Press). 


SER. 7, VOL. 43, NO. 345.—OCT. 1952 4A 


[ 1028 ] 


CII. The Decomposition of Thin Oriented Silver Halide Layers 


By D. W. PASHLEY 
Imperial College, London* 


[Received August 1, 1952] 


ABSTRACT 


Thin layers of oriented silver halides (AgCl, AgBr and AgT) have been 
prepared in various ways, and their stability to electron (30-60 kv) 
bombardment has been investigated. Unstable specimens are found to 
decompose to give highly oriented silver, in the case of silver chloride and 
silver bromide, but to give mainly randomly oriented silver in the case of 
silver iodide. Silver chloride also gives rise to highly oriented silver 
during ultra-violet irradiation. 

Chemical development of the halides produced randomly oriented silver, 
except when the halide was supported on a metallic silver crystal substrate. 
Tn this latter case the developed silver formed in an oriented fashion on 
the silver substrate. 

Diffuse streak patterns were observed from the silver halides and their 
origin is discussed. 


$1. INTRODUCTION 


Tue first definite evidence that the metallic silver lattice, formed during 
the decomposition of silver halide, can be oriented in a definite way 
relative to the halide lattice was obtained by Kooy and Burgers (1948). 
They studied, by means of x-rays, the silver formed during the photolysis 
of silver chloride, and obtained arc ring patterns indicating some preferred 
orientation of silver; this orientation was not identified. More recently, 
Berry (1949) and Berry and Griffith (1950) have studied in greater detail 
the effect of ultra-violet radiation on single crystals of silver bromide. 
They concluded that silver grows in a parallel orientation (two parallel 
cubic lattices) well inside the silver bromide, but that it has other orienta- 
tions at the surface. Simultaneously with this work, the decomposition 
of oriented silver halide layers during electron bombardment has been 
independently studied. Preliminary accounts of some of the results 
have already been given (Pashley 1950, 1951). The latter paper described 
how certain silver chloride specimens decompose to give metallic silver in 
an orientation parallel to that of the silver chloride ; prominent twinning 
of the silver on its (111) planes was also observed. 

The work has since been extended, and this paper describes a study of 
the decomposition of the three halides (silver chloride, silver bromide and 


* Communicated by Sir George Thomson, F.R.S. 
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silver iodide). Thin halide layers have been prepared in various ways, and 
their behaviour to electron bombardment examined by electron diffraction. 
It is found that some specimens decompose, and that others are quite 
stable. Some of the stable films have been decomposed by chemical 
means. A few specimens have also been reduced by ultra-violet 
irradiation. 


§2. EXPERIMENTAL PROCEDURE AND APPARATUS 


The silver halide layers were mostly prepared in one of the two following 
ways: (1) by evaporating the silver halide (in vacuo) on to various 
substrates, as described elsewhere (Pashley 1951, 1952 b) ; (2) by growing 
the halide by direct halogen attack of single crystals of silver (Pashley 
1952a). The substrates used for the evaporated layers included mica, 
rocksalt, potassium bromide, magnesium oxide and metallic silver, all in 
single crystal form. 

The specimens were always examined in the electron diffraction camera 
by the reflection technique ; by studying the diffraction patterns after 
various periods of electron bombardment (by the beam used to form the 
diffraction pattern), it could be ascertained whether the specimens were 
stable or unstable to the electrons. The unstable specimens were examined 
in detail, patterns being obtained and recorded at various stages of the 
bombardment, varying from about one minute to three hours. 

Some of the specimens were immersed in a dilute solution of a photo- 
graphic developer (metol-hydroquinone), and the growth of silver studied. 
No intermediate stages of growth could be recorded by this method, 
because the decomposition occurred too rapidly. 

The electron diffraction camera used was of a conventional type, 
employing a hot filament electron gun and magnetic focusing. The 
electron accelerating potentials were in the range 30-60 kv, and the camera 
length was 35 cm. 


§3. DECOMPOSITION OF CHEMICALLY GROWN LAYERS ON SILVER CRYSTALS 


Silver chloride, silver bromide and silver iodide films formed by direct 
halogen attack of single crystals of silver were completely stable to electron 
bombardment. These layers were usually highly oriented on the silver 
crystals (Pashley 1952 a), and have been decomposed by immersion in a 
chemical developer. It was found that the silver so produced was mostly 
oriented parallel to the substrate silver crystal, and was twinned on its 
(111) planes. There was also some randomly oriented silver. This type 
of silver formation was observed with layers of all three halides. The 
general features of the silver formation were independent of which lattice 
plane of the substrate silver was parallel to the surface (experiments were 
performed with (111), (100), (110) and (311) silver surfaces), and of the 
orientation of the halide layer. The extent of the twinning did vary from 
specimen to specimen. When the silver crystal had a surface containing 
many atomically flat patches, very prominent twinning occurred on (111) 
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planes parallel to the surface, but only relatively small amounts of silver 
were twinned on the other (111) type planes. In other cases, twinning 
was equally prominent on all four type (111) planes. It appears that the 
twinning occurred very readily on (111) planes which were parallel to 
facets bounding the substrate surface, but much less readily on other (1 11) 
type planes. Presumably the twinning does not occur until the particular 
(111) facets are developed. 

These effects show that the formation of silver occurs by deposition on to 
the substrate silver, either at the silver-silver halide interface, or via 
ionic conduction through the solution. This conclusion is supported by 
the fact that electrolytic deposition of silver on to the silver crystals gives 
rise to similar twinning effects. 


§4. DEcomMPOSITION OF EVAPORATED LAYERS OF SILVER BROMIDE AND 
SILVER CHLORIDE 


In order to eliminate the effect of the substrate silver, films of silver 
halides (thickness varying from 10 A to 1000 A) were prepared by evapor- 
ation on to cleavage surfaces of mica, rocksalt, potassium bromide and 
magnesium oxide. These layers were highly oriented (Pashley 1952 b), 
and were found to reduce to silver slowly during electron bombardment 
in the diffraction camera. Their decomposition could therefore be studied 
very conveniently without the need for any chemical developer. 


(a) Orientations during Electron Bombardment 


The electron beam was allowed to impinge on a certain spot 
of the specimen, and the growth of silver on this patch of the surface was 
studied as a function of time. The rate of decomposition varied markedly 
from specimen to specimen, and although thicker films decomposed more 
slowly than thinner films, the correlation between rate of decomposition 
and thickness was not as marked as was indicated by the first few experi- 
ments (Pashley 1950). 

The decomposition of silver chloride and of silver bromide occurred in 
very similar ways. A series of diffraction photographs showing the growth 
of silver during electron bombardment of silver chloride layers on mica 
has already been published (Pashley 1951). The main effects illustrated 
were that the silver chloride gradually decomposes to give a parallel 
orientation (two cubic lattices parallel) of silver, with prominent (111) 
twinning of the silver. Figure 1 (figs. 1-3, 6-8, Plate LX VIII) shows the 
pattern from a silver bromide layer on potassium bromide, after fifteen 
minutes bombardment. The smaller square pattern of very elongated spots 
is due to the silver bromide crystallites, which have fairly smooth tops. 
The larger square pattern of spots, which appears only after a few minutes 
bombardment, is due to the growth of a parallel orientation of silver, As 
the growth is continued, these latter spots split up into groups of four 
subsidiary spots, and the silver bromide pattern is reduced in intensity. 
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After ninety minutes bombardment, the silver bromide on the patch of 
specimen under bombardment had been almost completely decomposed, 
and the pattern of fig. 2 was obtained. The strongest pattern arises 
from the parallel orientation of silver, and forms a square pattern of spots. 
These spots are split into groups of four subsidiary spots; this splitting 
is explained in terms of the intersection of the Ewald sphere with the 
diffuse rods of intensity along all of the (111) type directions in the 
reciprocal lattice, and passing through all of the reciprocal lattice points. 
The pattern of subsidiary spots is identical with that obtained by Kirchner 
and Lassen (1935) from a layer of silver condensed on to rocksalt. These 
diffuse reciprocal lattice rods are due either to the presence of prominent 
(111) boundary planes, or to effects arising from repeated twinning on 


Fig. 4 
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(111) planes. Since the twinning probably only occurs after (111) facets are 
developed on the growing silver, it is likely that both effects (facets and 
twinning) are involved in the appearance of the subsidiary spots. The 
rings on fig. 2 are due to randomly oriented silver, which formed to a 
greater extent in silver bromide than in silver chloride. 
The occurrence of twinning was conclusively revealed when other 
azimuthal settings of the substrates were used. Figure 3 shows the pattern 
obtained from a silver bromide layer on rocksalt, after thirty minutes of 
electron bombardment. The interpretation of the pattern is given in 
fig. 4, there being five component patterns. The ( 111) twinning was a 
general characteristic of the silver growths which occurred during the 
decomposition of all of the specimens. In the very early stages of the 
growth the twin spots were either very weak, or absent altogether. 
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The subsidiary orientation which is listed in fig. 4 may be expressed as : 
(211) silver parallel to (100) silver halide, with [011] silver parallel to 
[011] silver halide, Again, this subsidiary orientation occurred for both 
silver bromide and silver chloride, when mica, rocksalt, and potassium 
bromide were used as substrates. It is therefore reasonable to assume that 
this orientation is not influenced by the substrate, but is a characteristic . 
of the decomposition process. Although the (211) silver planes and the 
(100) silver halide planes are parallel for this orientation, it does not 
necessarily follow that these are the contact planes between the two lattices. 
Table 1 lists the low-indexed halide planes which have a low-indexed 
parallel lattice plane in the silver. All these pairs of planes are possible 
contact planes. There are twelve orientations which are equivalent, as 


Table 1. Parallel Lattice Planes for the Subsidiary Orientation 


Halide Plane (100) (O11) (O1T) (111) (211). (211) (122) 
Parallel Silver Plane (211) (111) (O17) (011)_-—((122)_-—« (100)~—((21) 
Fig. 5 . 
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AgBr[i00} Ag[ori) eS / Projo 
ie 
AgBr [oil] - w x e 
Two rectangles represent two 
equivalent ways for Ag (011) 
to fit on AgBr (011). 


Fit of silver on silver bromide for subsidiary orientation, if (011) planes are 
contact planes. 


regards the fit between the two lattices, to this subsidiary orientation. 
For example, if the (011) planes are regarded as the contact planes, there 
are two equivalent ways in which the (011) silver plane can fit on the (O11) 
halide plane (see fig. 5). Since there are six non-parallel (011) type halide 
planes on which to fit, this makes twelve possibilities in all. Diffraction 
patterns from these equivalent orientations did not always occur when they 
would be expected, since a given one of the equivalent orientations only 
gave a strong pattern when the electron beam was set along a direction 
perpendicular to the parallel silver (111) and halide (011) type planes. 
At other settings, a given one of the equivalent orientations gave either no 
pattern (when one would be expected) or only a very weak pattern. This 
is not entirely understood, but possibly means that the silver which is in 
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the subsidiary orientation grows in the form of (11 1) plates (several layers 
thick) parallel to the silver halide (110) type planes, so that strong patterns 
will be obtained only when the plates present large areas to the electron 
beam. Thisisa possible explanation, since more electrons will impinge on 
a given plate when the electron beam is set perpendicular to it, than at any 
other setting. 

It is also interesting to note that no disorientation of the silver halide 
was observed during any stages of the decomposition. 

Several weak spots on the diffraction patterns were not interpreted. 


(b) The Influence of the Substrate 


In order to determine whether the substrate upon which the halide is 
deposited influences the occurrence of oriented silver, several different 
substrates have been used (see § 2). The results of the decomposition on 
all the substrates were entirely consistent, in that the observed orientations 
of silver were always the same relative to those of the parent silver halide 
(see § 4 (a)). In particular, the silver halide was in a ‘ fibrous ’ orientation 
on magnesium oxide, and the silver formed during decomposition was in 
the same ‘ fibrous’ orientation. 

A further experiment was carried out in order to obtain conclusive 
results on the influence of the substrate. Silver bromide could be oriented 
on potassium bromide in one of two ways, by changing the temperature of 
the substrate during deposition (Pashley 1952 b); (a) a parallel orient- 
ation ; (b) an orientation with (111) silver bromide planes parallel to the 
potassium bromide (100) surface. In both cases decomposition occurred 
such that the silver formed in the same orientations relative to those of 
the silver bromide parent crystal; the orientations of the silver relative 
to that of the potassium bromide substrate were therefore different in 
the two cases. There can thus be no doubt that the growth of silver 
described in § 4 (a) is characteristic of decomposition of the silver halide, 
and that the substrate does not initiate the orientation of the silver. 

Trillat (1951, 1952) has recently studied the decomposition (by electron 
bombardment) of oriented layers of silver bromide which were stripped 
from rocksalt substrates, and examined by the transmission method. 
He found non-oriented growth of silver. Similar transmission specimens 
have been examined by the author, with both silver bromide (stripped from 
potassium bromide substrates) and with silver chloride (stripped from 
rocksalt substrates). It was found that although much of the silver 
formed during bombardment was randomly oriented, there was also some 
parallel oriented silver, particularly in the case of silver chloride. Slight 
parallel orientation of silver has since been obtained by Trillat with his 
silver bromide specimens (private communication). This result adds 
further confirmation to the conclusion that the parallel orientation of 
silver is a characteristic of the decomposition process, and is not just a’ 
substrate effect. A further important effect (discussed in §8) was that 
the silver halide recrystallized during electron bombardment, so that some 
of it became randomly oriented before being decomposed. 
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(c) Decomposition by Ultra-violet I rradiation 

Preliminary experiments have been carried out on the decomposition 
of evaporated silver halide layers during ultra-violet irradiation. Silver 
chloride was evaporated on to some mica and rocksalt specimens, and two 
of the specimens were used to check that the chloride had oriented on the 
substrates in the usual manner. The other specimens were then placed in 
an evacuated quartz tube, and irradiated by means of an ultra-violet 
lamp for eight hours. These specimens were not subjected to electron 
bombardment before the irradiation. After irradiation, the specimens 
were examined in the electron diffraction camera ; the pattern obtained 
immediately the electron beam was switched on indicated that much 
decomposition had occurred during the ultra-violet irradiation. The 
pattern from such a decomposed layer on rocksalt is shown in fig. 6. 
The growth is identical to that which occurred during electron bombard- 
ment (the larger square pattern of fig. 6 indicating the parallel orientation 
of silver), except that rather more randomly oriented silver is formed than 
in the electron case. The study of the layers on mica confirmed that the 
decomposition of silver chloride under ultra-violet irradiation occurs in a 
manner which is not significantly different from that occurring during 
electron bombardment. 

Some further experiments, in particular with silver bromide, have 
indicated that more complicated effects can occur, but since some of these 
effects might be influenced by reactions with the substrate, they will not 
be reported at present. Itis worthy of note, however, that in none of these 
experiments have the non-parallel orientations of silver reported by Berry 
and Griffith (1950) been observed. 


(d) Chemical Development 


Specimens prepared on magnesium oxide substrates have been 
decomposed by treatment with a dilute chemical developer (metol- 
hydroquinone), without subjecting the specimens to any previous electron 
bombardment. Subsequent electron diffraction examination of these 


specimens has in all cases revealed the formation of randomly oriented 
silver. 


§ 5. Decomposition oF EvaporatED LAYERS OF SILVER IODIDE 

Silver iodide layers have also been prepared by evaporation on to various 
substrates. Although the molten bead of silver iodide, from which the 
evaporation was effected, gradually decomposed, the diffraction patterns 
from the freshly prepared specimens were due entirely to silver iodide. 
There are two common forms of silver iodide, a cubic (zinc-blende type) 
structure and a hexagonal (zincite type) structure. Both modifications 
often occur together. It was found, however, that a layer consisting of 
only one modification could be formed. Such layers were highly oriented. 

Figure 7 shows a pattern from an oriented layer of hexagonal silver 
iodide on mica, after twenty minutes exposure to the electron beam. Most 
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of the silver formed has grown in a random orientation, giving rise to the 
ring pattern. The spot pattern arises from the oriented silver iodide. 
With continued bombardment, the silver pattern becomes stronger, and 
many weak arcs can be distinguished on the silver rings. Two distinct 
patterns have been identified, indicating the orientations given in table 2. 


Table 2. Silver Orientations in Silver Todide 
: Orientation (1) (2) 
Silver plane parallel to silver iodide (110) (110) 


Hexagonal | _(001). 
Iodide. |Silver axis parallel to silver iodide [110] [112 


[100]. 
Silver plane parallel to silver iodide (110) (110) 
Cubic mel): 
Todide. Silver axis parallel to silver iodide [110] (ti 
[110]. 


We can therefore say that the decomposition of silver iodide, both cubic 
and hexagonal, leads mainly to randomly oriented silver. It is not yet 
certain whether the weak orientations observed are influenced by the 
substrate, or not, but it seems probable that they are not. 


§6. THe STapiity oF SinvER HaLipEs TO ELECTRON BOMBARDMENT 


A point of interest that arises from the present work concerns the stability 
of silver halides to electron bombardment. The only specimens which are 
found to be unstable are: (1) evaporated layers of halide on rocksalt, 
potassium bromide, magnesium oxide and mica substrates ; (2) stripped 
layers from rocksalt and potassium bromide substrates. The following 
specimens showed no indication of decomposition after two hours’ bom- 
bardment in the diffraction camera: (1) chemically grown layers on silver ; 
(2) evaporated layers of halide on metallic silver. Layers formed by 
evaporating silver halide on to a formvar film were found to decompose 
slightly during a two-hour bombardment (by transmission). Several 
other workers have found that various kinds of silver halide specimens are 
completely stable to electron bombardment. A summary of these results 
has already been given (Pashley 1951). 


§7. DirrusE BAcKkGROUND PATTERNS FROM THE SILVER HALIDES 


Burgers and Tan Koen Hiok (1946) have studied silver chloride 
crystals by x-ray diffraction, and have observed the presence of diffuse 
streaks on their diffraction patterns. Similar streaks are also found 
on the electron diffraction patterns from some of the silver chloride 
and silver bromide specimens used in the present work. The streaks 
were more prominent with silver bromide specimens than with those 
of silver chloride. The streak patterns correspond to the presence of 
diffuse rods of intensity passing through all the ‘strong’ reciprocal 
lattice points (allowed points with even indices) along the three cube 
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directions. The occurrence of the streaks on patterns from sateen, 
grown layers on silver crystals is well illustrated on a photograph aaer < 
previously (Pashley 1952 a, fig. 13). The streaks from an evaporate 

layer are shown in fig. 8. None of the diffuse reciprocal lattice rods 
lies in the (211) reciprocal lattice plane, which is perpendicular to the 
electron beam, but some are inclined to it at only a small angle. The 
disposition of the diffuse rods just in front of the (211) reciprocal lattice 
plane (a similar set is also behind) is indicated in fig. 9. There are two 
orientations of the bromide, each of which gives the same spot pattern; the 


Fig. 9 


{010} ae on 


Cube edge directions for one 
silver bromide orientation 


fool} 
Cube edge directions for | 
other silver bromide orientatio 


‘Strong’ allowed spots in ‘Weak’ allowed spots in 
@ (211) reciprocal lattice plane. (211) reciprocal lattice plane. 
Disposition of diffuse rods in front of (211) reciprocal lattice plane, for two 
orientations. 


diffuse streak patterns from the two orientations are not coincident, the 
two components being indicated in figs. 9(@) and 9(b). When the electron 
beam is set nearly perpendicular to the (211) reciprocal lattice plane, 
the Ewald sphere intersects these broad diffuse rods to give the pattern 
plotted in fig. 10. This is identical with the diffuse pattern in fig. 8. 
Such diffuse streaks were observed quite frequently, but varied in 
intensity from specimen to specimen. It is considered that the reason 
for this variation is due to the variation in general background intensity 
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from specimen to specimen, since even quite a low background is 
sufficient to mask the diffuse streaks. 

Similar diffuse streaks were also found with evaporated thallium 
chloride layers (Pashley 1952 b), and can also be seen on an electron 
diffraction photograph published by Yamaguti (1935). These latter 
streaks were obtained from a cleavage surface of zine-blende from which 
a previously grown oxide layer had been removed. Yamaguti did 
not comment on the streaks; these streaks can also be explained in 
terms of the presence of broad diffuse rods of intensity along the cube 
directions of the zine blende reciprocal lattice. 

Burgers and Tan Koen Hiok (1946) suggested that the streaks which 
they observed could be explained by supposing that various cube planes 
of their crystals had slipped, so that they were not phase linked with the 
rest of the crystal. Mitchell (1949) has suggested that. the presence 


Fig. 10 


e 8 o 8 ° 8 ° 8 e 
egies Ser 8 se 3. 6 
Diffuse streak pattern of fig. 8. 
of F-centres could explain the streaks. The streaks observed on the 
electron diffraction patterns from the thin silver halide layers cannot 
be due to F-centre aggregates. A (100) plane of the halide, with its 
halide ions replaced by F-centres over a small area, would give rise to 
diffuse rods of intensity along [100] directions in the reciprocal lattice, 
and these rods would pass through all allowed halide reciprocal lattice 
points (points with indices all even or all odd). ‘Weak’ reciprocal 
lattice points (allowed points with odd indices) should therefore have 
rods passing through them of the same intensity as those passing through 
the: ‘ strong’ points. Streaks through the ‘ weak’ points were never 
observed. Since Burgers and Tan Koen Hiok found that their diffuse 
x-ray patterns were present even with the chloride at liquid air 
temperature, it seems unlikely that the streaks are of thermal origin. 
It is thus considered that the possible explanations of the streaks 
are: (1) the slipping of various cube planes, as suggested by Burgers 
and Tan Koen Hiok ; (2) the presence of numerous small cavities in the 


crystals, consisting of missing cube planes (both silver and halide ions) 
of small area, and one layer in thickness ; (3) some distortion involving 


the cube planes. 
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The presence, or otherwise, of the diffuse streaks on the electron 
diffraction patterns showed no correlation with the stability of the 
halide to electron bombardment. The streaks do not, therefore, indicate 
the early stages of the silver growth, or the property of a particular 
specimen which controls its stability to electrons. It is possible, however, 
that the property which gives rise to the streaks does play some vital 
part in the decomposition process. 


§8. SUMMARY AND DISCUSSION 


The main results presented above may be summarized as follows : 

(1) The decomposition, during electron bombardment, of thin layers 
of oriented silver bromide and silver chloride gives rise to a parallel 
orientation of silver, with prominent (111) twinning. Other weak 
orientations also occur. 

(2) Silver chloride behaves in the same way during ultra-violet 
irradiation. 

(3) Silver iodide decomposes during electron bombardment to give 
mainly randomly oriented silver. 

(4) Diffuse streak patterns are obtained from silver bromide and silver 
chloride; they do not appear to be linked with the early stages of 
nucleation of silver. 

(5) Chemical development of oriented layers of silver bromide and 
silver chloride gives randomly oriented silver except when a single 
crystal of silver is used as substrate; then the silver which is formed 
by development continues the structure of the substrate, and twins. 
on (111) planes. 

The orientation of silver produced during electron bombardment 
of silver halides has been shown (see § 4 (b)) to be a characteristic of 
the decomposition process, and is not a substrate effect. The apparent 
discrepancy between the results of transmission experiments (which 
indicate only slight orientation of silver) and the reflection experiments. 
is considered to be a due to a combination of the following: (1) the 
silver halide crystallites on the substrates are rigidly held, whereas the 
silver halide crystallites in the stripped films are not ; (2) recrystallization 
of the halide occurs with the stripped films (giving randomly oriented 
silver halide), but no disorientation occurs with the halides on the 
substrates. This recrystallization is possibly due to thermal effects. 
of the electron beam, for although transmission specimens might be 
appreciably heated (cooling conditions are poor), crystallites on a massive 
substrate will not be appreciably heated because of the conduction 
losses via the substrate. Thus the role of the substrate in the above 
experiments is to orient the silver halide crystallites rigidly, and to provide 
a means of keeping them cool, without appreciably influencing the 
orientation of silver produced during decomposition. 

Berry and Griffith (1950) find that, during the ultra-violet irradiation 
of a single of silver bromide, parallel oriented silver grows inside 
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the crystal, and that other orientations occur at the surface. The 
evaporated layers of halide which have been used in the present work 
have a very high surface to volume ratio, and it would be expected 
that orientations characteristic of surface growth would predominate. 
Instead, the observed orientations of silver in these layers agree with 
the orientations which Berry and Griffith find inside their crystals. 
The relation between these two sets of results is not clear. 

It is of interest to consider how the above results fit in with 
mechanisms suggested for photographic latent image formation. Although 
the original growth centres of the silver have not been detected by 
electron diffraction, the fact that highly oriented growth of silver occurs 
during electron bombardment of silver chloride and silver bromide 
shows that the growth centres themselves must be well oriented inside 
the silver halide lattice. In this connection, a latent image theory 
proposed by Huggins (1943) is of interest. Huggins suggests that the 
normal rocksalt structure of silver bromide might be converted to a 
zinc-blende type of structure when photoelectrons are released inside 
the crystal. The zinc-blende structure is then supposed to break down 
to form metallic silver. If such a mechanism were to occur, then electron 
diffraction would be expected to reveal the presence of the abnormal 
structure. No evidence of this structure was ever observed. 

The mechanism by which the decomposition takes place must be such 
that : (1) oriented silver nuclei are formed ; (2) these nuclei can continue 
to grow, with the orientation of the silver and the unchanged halide 
maintained; (3) the liberated halogen atoms can freely escape; 
(4) twinning of the silver can occur ; (5) almost complete decomposition 
can be achieved, without destroying the orientation. It is difficult 
to see how such conditions of growth can be maintained unless either 
(a) the reaction takes place at the substrate-halide interface, or 
(b) a certain amount of silver halide remains unchanged, and acts as 
a support for the oriented silver. Since it has been shown that the 
substrate is not responsible for the initiation of the orientation, (b) seems 
the more likely. Further, there were always at least small traces of 
silver halide indicated on the diffraction photographs, and in most 
cases (except for very thin layers) unchanged silver halide very close 
to the interface with the substrate would not be detected in the electron 
diffraction pattern. 

Whether the nucleation of the silver occurs at the surface, or inside 
the crystallites, cannot be determined directly. If, however, nucleation 
started at the surface, it might be expected that the orientation of 
the silver would be dependent on the crystal facet upon which nucleation 
took place, just as in the case for the growth of silver halides on silver 
(Pashley 1952 a). No variation in orientation from specimen to specimen 
was found, despite the fact that the halide layers on rocksalt and mica 
almost certainly had different facet structures. It therefore seems more 
likely that the nucleation was not a predominantly surface phenomenon, 
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and that much of the growth started inside the halide crystallites. This 
is, at least, in partial agreement with the results of Berry and Griffith 
(1950), who find that the parallel growth occurs inside a massive silver 
bromide crystal. 

It would be of interest if any distinction could be made between 
the original Gurney and Mott (1938) theory of latent image (involving 
growth of silver on to ‘sensitivity specks’, via motion of interstitial 
silver ions) and the modification of this theory proposed by Mitchell 
(1949). Mitchell suggests that a latent image in a silver bromide grain 
could be an aggregate of F-centres ; this aggregate can grow to a certain 
size, and then collapses to form a metallic silver nucleus. A mechanism 
of this kind might well be expected to form a parallel orientation of 
silver, since the face-centred silver ion lattice in the halide structure 
would then have to contract, without change of orientation, to form 
the face-centred cubic lattice of metallic silver. However, this is not 
a unique interpretation, and the occurrence of twinned silver definitely 
indicates that motion of silver ions, either as a surface mobility or as 
a mobility of interstitial ions, does occur after a certain stage of the 
process. 

It is thus clear that there is not sufficient evidence available at present 
to decide on the detailed mechanism whereby oriented silver occurs 
in silver chloride and silver bromide (with hardly any orientation in 
silver iodide). However, the experiments described here have given 
a set of results which must be explained by any comprehensive theory 
of latent image formation. 
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EXPLANATION OF THE PLATE 


Fig. 1 
Silver bromide layer on potassium bromide, after 15 minutes electron 
bombardment. Bromide [001] azimuth of (100) surface. 


Hig. 2 
As for fig. 1, except 90 minutes electron bombardment. 
Fig. 3 
Silver bromide layer on rocksalt, after 30 minutes electron bombardment. 
Bromide [011] azimuth of (100) surface. 
Fig. 6° 
Silver chloride layer on rocksalt, after ultra-violet irradiation for 8 hours. 
Chloride [001] azimuth of (100) surface. 
Fig. 7 
Hexagonal silver iodide layer on mica, after 20 minutes electron bombard- 
ment. Iodide [100] azimuth of (001) surface. 
Fig. 8 


Silver bromide layer on mica, after 5 minutes electron bombardment. 
Bromide [211] azimuth of (111) surface. 


[ 1042 ] 
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ABSTRACT 
A discussion of the virial theorem for both Thomas—Fermi and Thomas— 
Fermi—Dirac systems is given, the results being more general than have 
hitherto been obtained. 
Finally the particular case of a homonuclear diatomic molecule is 
briefly considered. 


$1. INTRODUCTION 


THE virial theorem, as is well known, holds in classical and quantum 
mechanics and in both theories it has had useful applications. 

However it is by no means obvious that the theorem will remain valid 
in the Thomas—Fermi (T.F.) theory in view of the simplifying assumptions 
involved. Fock (1932) made the first investigation on the subject. 
He considered a T.F. system of many electrons and nuclei, and using a 
variational method he derived the form of the virial theorem for such 
systems. His proof was however later shown by Jensen (1933) to be valid 
only for the case of a single centre, on account of the improper variations 
employed. Jensen attempted to avoid the difficulties of Fock’s treatment 
but in doing so made some restrictive assumptions which subsequently 
proved to be unjustified. Duffin (1935) gave the first satisfactory 
derivation of the theorem for T.F. systems of many electrons and nuclei. 

The proofs mentioned previously all deal with the original T.F. theory 
in which exchange is neglected. However Jensen (1934) has shown how 
the theorem is modified by the exchange term for the case of a single 
nucleus, using Fock’s method. 

As was first pointed out by Slater and Krutter (1935) the T.F. method is 
of considerable value for discussion of the behaviour of material under 
high pressure. Hence it is of interest to investigate the form of the virial 
theorem when the systems considered are not free, but have a pressure on 
the boundary. For the case of a single nucleus such an investigation has 
been carried out by Slater and Krutter (1935) and Jensen (1939). 

Recently, an ingenious discussion of the virial theorem based on 
similarity considerations has been given by Feynman, Metropolis and 
Teller (1949). 

However, no unified and general proof of the theorem has yet been 
given, and indeed in the Thomas—Fermi-—Dirac case the form of the theorem 
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in the many nuclei problem has not been established. The purpose of 
this work is to supply such a general treatment from which all existing 
results can be obtained immediately as particular cases. 


§2. THe Unmopiriep THomas—Frermt Mernop 

The treatment that we shall adopt follows that given by Duffin (1935), 
and hence we shall outline his method here. We shall see that we can 
immediately obtain the required generalization to systems which are not 
free, simply by retaining a surface integral which Duffin quite justifiably 
dropped for the free systems that he was considering. 

We take as starting point in all cases the expression for the kinetic 
energy 7’. In the original T.F. method this may be written in terms of the 
total electrostatic potential V as 


3 
Tee eer beds wareee Atay 
32m 
where ee oe (2me)?”° 


and —eV, is the maximum energy of an electron. 
We now use the relation , 


3 | (V—V,)@dr= | We er. ds | (VV er ata Vda. (2) 
os s a 

Duffin now drops the first term since for free systems V—V,=0 on the 

boundary. This is not true however for systems in which there is a 

pressure on the boundary, and hence we retain the term here. Using the 

usual T.F. expression for the charge density p in (2) and then substituting 

in (1) we obtain 


a = ih |. (VV, r-ds-+ | pr. grad Vd BOSD) 


Ife, and r,, are the charge and position vector of the kth nucleus, and V, is 
the potential due to the electrons, then we can write 


(pe eee eee ml (Y) 


k |r—r,,| 


Also Duffin has proved the following lemma : 


Jor-srad V.dr=—3| p Vi Gigi ee ae 2. OD) 
and hence 
Ei 5/2 Ck 
oT = fe. (V=V, 2 r. ds—3| pV.dr+2| p regrad (“tS de (6) 
We now use the identity 
r. grad [r—r,,|-!=—|r—r,|[-!—r,- grad, |r—r,[7 - . - (7) 


and proceeding exactly as Duffin has done we arrive at the result 
ee ok LS) yp sds lim Séethne e; (te) s eee 8 
eg | iro rods him 4 eg ry €q (Fx) (8) 
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where U is the potential energy of the electrons and «, is the electric field 
produced by all charges outside small spheres of radius a surrounding the 


unclei. 
Evidently this is the general form of the theorem, the first term on the 


right arising from the fact that the system considered is not free, in other 
words it is to be interpreted as the virial of the forces acting on the 
boundary, whilst the second term can be interpreted as arising from the 


forces required to hold fast the nuclei. 
The various forms of the theorem obtained hitherto follow immediately 


from (8). Duffin’s form is obtained as we have said, by putting the first 
term on the right equal to zero. Further, for the case of one centre the 


second term on the right is obviously zero. 
For free atoms and ions (V—V,.)=0 


on the boundary, and hence 2T+U=0. 
‘For non-free atoms we have 


2T+U = a (V—V,)?/* r. ds 


= oZ ( Vio). femenaare . 30, 


where v is the volume. 
But Slater and Krutter and Jensen have shown that the pressure p is 
simply given by 


at as . 
1 10n [Ch == Vig) herendace 
and hence 2T+U=3pv. 


§3. THe THomas—FrermiI—Drrac Mrernop 


The preceding proof will now be generalized to apply to the Thomas— 
Fermi—Dirac (T.F.D.) method in which exchange is included. The only 
discussions given hitherto (Jensen 1934, Feynman, Metropolis and Teller 
1949) are limited to the case of a single centre. 

As before, we take as starting-point the expression for the kinetic energy 
fe 


T=c, | n! dr, Po, a. ae a 
where = oe oe 
“k= 10m \8nr 


and 2 is the number of electrons per unit volume. For T.F.D. systems, 2 
is related to the total electrostatic potential V by 


n= 2 [a+ (V—Vota)irp, eee MeN ee 


(2me3)1/2 


h = 
where a ; 
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It is convenient to introduce the abbreviation 


Ue (VV te yep e ..  e) 
and then from (9) 
Su 
Das | 058 ar. 55s ly SARS Re as eR lb} 
Using Gauss’ theorem this can be written 
" ¢ i 5/2 
or = | Ys r.ds+ | pr. grad U dr. Stee ae 10'S) 
Now grad U=grad V+2a grad U!/2 see sit, se 14) 


and hence 
27 = i | Ue r. ds+| p r.grad Vdr— |v r.grad U1? dr. 
aT 2a 
The last term can evidently be written as 
= ae . grad U? dr, 
and using Gauss’ theorem this becomes 


—E([o- r.ds—3 | 0 ir). 


The term involving the volume integral can be shown immediately to 
equal the negative of the exchange energy A, 


_ where A=—c, | nt dr 
e2 31/3 
and a eT Te 


Proceeding entirely as before we obtain 
é 5a 
f+U+A= | 32 — —_**)r.ds+Lim +e, r,.«, (r,) - . (15) 
107 4 a=0 k 
or writing the first integral in terms of 


2T+U-+A = | (Foun 9—4e,n8?) Kas onde NOE r,-€, (;,)- « (16) 


This is evidently then the general form of the virial theorem for T.F.D. 
systems. 

As before we can immediately obtain all previous results as special cases. 
For the case of a single centre the second term on the right is zero. Also 
Jensen has shown that in this case, the pressure p defined by 


a GE 
P=— dy 
is given by p=(3e,n?— 3c ,n4!),oundary 
and hence 2T +U+A=3pv. SOE PT Ss Sree ame Ld) 


This is the result first obtained by Jensen. 
4B2 
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§4. Case or A Frez HoMONUCLEAR DIATomIc MOLECULE 
The results obtained previously are of course quite general, applying to 
any molecule or crystal. In concluding the discussion, we wish to point 
out that the theorem in the special case of a homonuclear diatomic 
molecule described by the original T.F. theory can be put into the more 
usual form* 


dE 


27-+U=—R =, (18) 


where R is the internuclear distance. 
This result, first given by Slater (1933), has not hitherto been proved to 
hold in the T.F. theory as far as the writer is aware. 
All that is necessary is to discuss the last term in (6), that is 
— 5 | ne (« . grad we) dr. 
k [r—r;.| 
In the case we are considering this can be written in the form 


| meVxedr+ | ner oe dr, 


where Vy is the potential due to the nuclei. We have now to set up the 
energy expression and differentiate it with respect to & to obtain 
dE dV x 
dR =—eé | n ahi dr. 
Hence eqn. (18) is established. It can be shown, after some 


manipulation, that in this case the kinetic energy 7 and the potential energy 


U can be written in the form 
6 _ dH 
Pig erRe 
12a 5 dE 
U= 7 Ze(h dn tz hap, 


T=— “Ze V dy 
(19) 


where (V,)y is the potential at either nucleus due to the electrons. 
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ABSTRACT 


The volume of vapour evolved from liquid helium by known amounts 
of heat has been measured at temperatures between 2-9 and 4:5°xK. 
The derivation of latent heat values from these results is discussed, and 
it is shown that the existing experimentally determined values of the 
saturated vapour density are not sufficiently reliable for this purpose. 
The probable true values of the vapour density are considered. 


§1. INTRODUCTION 
WE have constructed an apparatus for measuring very small heat 
capacities and small differences of heat capacity. In order to avoid 
errors caused by the heat capacity of the calorimeter in the conventional 
type of apparatus we use an improved version of Dewar’s apparatus (1905) 
in which the heat content of the sample is measured by determining the 
volume of vapour driven off by this heat from a low boiling liquid. 
From this a mean specific heat can be derived. 

In our experiments we measure the increase in volume of the liquid— 
vapour system at constant pressure ; but this increase in volume repre- 
sents a smaller mass of gas than is produced from the liquid. Suppose 
1 g of the liquid to be evaporated ; the gramme of vapour produced can 
‘be considered in two parts: (1) the vapour which replaces the liquid 
evaporated, and is of mass py/p; g and (2) the vapour of which the volume 
at room temperature is actually measured and is of mass (l—py/pz) g. 
pz, and py are the liquid and vapour densities. 

In order to calculate the number of moles of liquid evaporated we 
must divide the measured volume of gas by (l—py/p,). In the case of 
helium above 3° kK this factor is appreciably less than unity (~0-86 at 
4-2° x) and cannot be calculated accurately owing to uncertainty as to the 
density of the saturated vapour. This uncertainty is due to the great 
deviation from ideality of the vapour and there are few direct measure- 
ments of the density at saturation. 

As the only measurements of the latent heat of helium at different 
temperatures are those of Dana and Kamerlingh Onnes (1926), in which 
the volume of liquid evaporated was determined, we have measured 
directly the apparent latent heat, L,, in terms of gas evolved, which is 
the quantity actually required for our specific heat experiments. This 
is related to the true latent heat, L, by the equation L=L,(1—py/p;). 
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§ 2. EXPERIMENTAL 


1. Apparatus 
; re used was the one designed primarily for the heat or 
measurements. The relevant part of the apparatus is shown in fig. " 
A narrow (2mm diameter) thin-walled, German silver tube fe aes ‘ 
the annular copper vessel A to a copper tube passing throug Be . nig 
pressure pot 8, of a Simon-type helium liquefier, and then ) oa 
source of helium gas and to an apparatus for making continuous measure- 
ments of volume changes. 
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Apparatus. 


The pressure in the system is measured on a mercury manometer, and 
by trapping this pressure in a litre flask C, the differential oil manometer D 
can be used as a sensitive detector of pressure change. The flask C is 
immersed in water to keep its temperature constant during the course of a 
single experiment (about 30 minutes). E consists of a tube of accurately 
constant cross-section, and is actually one of three ‘ Veridia’ tubes of 
different diameters, used to allow various volume ranges to be measured. 
By changing the level of the mercury in K, the volume of the system can 
be changed while keeping the pressure accurately constant. 
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2.2. Method 


Helium is liquefied in S, and then helium gas let into tube B to allow 
it to be condensed into A. The level of the mercury in E is pushed by 
compressed air up to the highest possible level. 

On account of stray heat influx a steady flow of vapour leaves A and is 
measured by lowering the level E at such a rate that the pressure in the 
system is constant ; the height being read on a glass mirror scale. The 
height of mercury is read off every minute for 6 to 10 minutes. The points 
are plotted as soon as they are taken so that we can judge the degree of 
linearity of the background at any time. 

Heat is then supplied to A by means of an electric heater H, for a 
measured time. On switching off the current the background is again 
reached after a few minutes. The vertical ‘ height’ between initial and 
final backgrounds then gives the volume corresponding to the heat 
supplied by the heater. 

Wexler (1951) has carried out a similar determination at 4-228° K. 
In his experiment he found the volume evolved to depend upon the rate of 
heating as well as on the quantity of heat supplied. The large volume of 
escaping cold vapour reduced the temperature gradient in the equivalent 
of tube B and thus reduced the appreciable background caused by heat 
conducted down this tube. In our apparatus the background is almost 
entirely due.to radiation down the wide central tube F, provided as a 
means of dropping specimens into the helium region, and not to conduction 
along tube B. To verify this, we measured the background evaporation 
of helium when the excess in pressure on S over that on A was increased 
from the usual 1cm of mercury to 15cm. The background changed 
from 0-2 cm? of gas per minute to 0-4 cm?, showing that under our normal 
conditions only about 1/15 of the total background is due to conduction 
down tubes B and F. Any loss in background due to a reduced temper- 
ature gradient in tube B should then, for the times and rates of heating we 
used (generally ~30 cm? gas evolved in ~2 minutes), produce an error in 
the latent heat of less than 1/10°%. We did indeed find that the values of 
latent heat were independent of the rate of heating. 

Determinations of latent heat at pressures lower than normal atmo- 
spheric pressure are carried out by pumping the helium in §, and allowing 
the pressure in the main system to follow. As the pressure in 8 must be 
slightly in excess of that in A, to enable evaporating vapour to escape, 
once the pressure in A has reached the required value the rate of pumping 
is decreased and adjusted to stabilize the pressure in 8 at the required 
value. When higher pressures are required, the gas evaporating from 
Sis made to pass through one of a series of fine capillaries chosen to give 
the required pressure heads. In this way constant pressure heads up to 
about 20cm of mercury were applied. The pressure in the system is 
allowed to build up to the required value, and then the experiment is 
carried out in the same manner as at lower pressures. 
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§3. Resutts AnD DiscussIoNn 


3.1. Experimental Results 

The ‘ apparent ’ Jatent heat of helium was measured between 2-9° and 
45°. The values found are tabulated in units of calories per mole of 
gas collected (table 1). At 423° we find agreement with the value 
determined by Wexler. 

We calculated the true latent heat at 4-23°k from our values of the 
apparent latent heat by using the values of the liquid density measured by 
Kamerlingh Onnes and Boks (1924) and of the saturated vapour density 


Table 1. Experimental Values of Latent Heat 


Apparent True 
Temperature | latent heat ears latent heat 
“kK L£,+0-1 ashe L,(1—py/pt) 

cal/mole cal/mole 
2-867 22-9° 0-976 22-4 
3:172 23-4 0-958 22-4 
3-180 23°3 _ 0-958 22-35 
3-714 235 0-922 21-6 
3°987 23°25 0-893 20-75 
3-999 23°3 0-892 20°8 
4-198 22-80* 0-865 19-72 

+0-05 

4-280 22-6 — = 
4-280 22-65 oo — 
4-434 22-2 —s a 
4-440 22-1 ats ma 
4-510 21-6 = = 
4-512 PA ibe — eo 


Extrapolated values at 
normal boiling point 


* Mean of six values. 


All calculations are made using J =4:18 joules/calorie. 


measured by Mathias, Crommelin, Kamerlingh Onnes and Swallow (1925) 
both taken from Keesom (1942). As was found by Wexler, this calcu- 
lated latent heat does not agree with the direct determination by Dana 
and Kamerlingh Onnes and is about 1% too high. However, at the 
lowest temperatures at which we have made measurements ihe factor 
(1— pv/px) is near unity and can be calculated with more certainty from 
the virial coefficients. At these temperatures, by calculating the vapour 
density from the first three virial coefficients, using the table of ‘adopted ’ 
values given by Keesom (1942), we find good agreement between our 
corrected latent heat and the values of Dana and Kamerlingh Onnes. 
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This suggests that the discrepancy at higher temperatures is due to an 
error in the calculated correction factor, and therefore probably in the 
measured vapour density. 

In order to obtain other values of the saturated vapour density we have 
plotted the isotherms, pressure versus density, of helium vapour obtained 
by Keesom and Walstra (1940). These extend to about 90°% of the 
saturation pressure and we have extrapolated these up to the saturation 
pressure and read off the density. It was found that the experimentally 
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4 Calculated from virial coefficients. 
I. Saturated vapour density of helium +} @ From extrapolation of isotherms. 
O Experimental (Mathias e¢ al.). 


II. Vapour density which helium would have if it were ideal at saturation 
pressure. 


determined saturated vapour density at 4-23° K is incompatible with the 
extrapolated isotherms. We have calculated some additional values of 
the saturated vapour density below 4° K by means of the virial coefficients 
given by Keesom and Walstra. Figure 2 shows the vapour density 
derived by these two methods, and also the direct measurements of 
Mathias ef al. 
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Using values of vapour density taken from the curve of fig. 2 we have 
calculated the true latent heat from our experiments, as shown in table 1. 
We obtain good agreement with the measurements of Dana and Kamer- 
lingh Onnes as can be seen in fig. 3. Our results go to higher temperatures 
than do those of Dana and Kamerlingh Onnes, but the uncertainty m 
py prevents us from converting our measured values into true latent heats. 
In this region any error in py would cause a serious error in the calculation. 


Cals/mole 


2S SI@) eis 4-0 4°5. 
Temperature °K 


; — Dana and Kamerlingh Onnes. 
Eaton aot { ® Present experiments. 


3.2. Comparison with Calculated Latent Heat 


Further evidence that the saturated vapour density is in error is 
provided by the discrepancy between the measured latent heat and the 
value calculated by Kistemaker (1946) from the Clausius—Clapeyron 
equation, dP/dT’=L/T(1/py—1/p,). In this calculation any error in the 
vapour density has a greater effect on the value of the latent heat derived 
than it has on the value derived from our measurements, and this is 
indeed found. 

We have calculated the latent heat for various temperatures from the 
Clausius—Clapeyron equation using the values of vapour density from 
fig. 2. Table 2 shows the data used and the results of these calculations - 
it is seen that there is agreement with the experimental values within the 
accuracy obtainable from the data used. The table includes values for 
one temperature as given by Kistemaker. It will be noted that some 
discrepancies occur due to the different vapour pressure—temperature 
scales used. We have taken values from the Royal Society Mond 
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Laboratory’s tables (1949), whereas Kistemaker uses the ‘1937’ scale 
to which he has applied corrections calculated from his own vapour 
pressure experiments. 

Below 4-2°k Kistemaker has calculated the latent heat by equating 
the Gibbs free energies, G, of the liquid and of the saturated vapour for a 
particular temperature. At 3-5°xK and below, these calculated values 
agree with the experimental values, but at 4° and 4-2° x the discrepancy 


Table 2. Calculation of Latent Heat from Clausius-Clapeyron Equation 


dP /1 1 
Temperature dP/dT x 106 <1 oe e(en63 | fh ( ) L exp. 


eens dT \ py 
: PV pL | 
K cg units+2 0-01 / allaole cal/mole 


3°10 395 


0- 

Sane, 676 1: 
4-23 979 1: 

Kistemaker’|s calculation] :— 
4-23 996 


* See footnote. 


is of the order of 4 to 5%,7 whereas the stated uncertainty of the 
calculations is only 14%. However, the virial coefficients which he uses 
at 4° kK give a vapour density approximately in agreement with fig. 2. 
Thus the discrepancy must be elsewhere in Kistemaker’s calculation by 
this method. 


§4. SuMMARY 


We have determined the latent heat of helium between 2-9° and 4:5° k 
by measuring the volume of vapour produced by known amounts of heat. 
Using values of the saturated vapour density found by extrapolating 
isotherms we obtain good agreement with the experimental values of 
Dana and Kamerlingh Onnes between 2-9 and 4-2°K. _ This suggests that 
the discrepancies found by Kistemaker between the experimentally 
determined and the calculated values for latent heat at 4° and 4-2° kK are 
not due to an error in the experimental values, as Kistemaker thought 
probable. 

We have not converted our measurements above the normal boiling 
point into true latent heat as we do not think that the experimentally 
determined saturated vapour densities are sufficiently reliable. 

+ It appears that the ‘smoothed’ value of the experimental latent heat 
at 4-25°K given by Keesom (4.58 cal/g) is a misprint and should actually be 
4-85 cal/g, and also that the value at 4:2°K given by Kistemaker is in error 
due to this misprint. 
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Values of the latent heat calculated from the Clausius—Clapeyron 
equation agree with the experimental values, but again the accuracy is 
limited by the lack of knowledge of vapour densities. 

In view of the revealed uncertainty in vapour density values it is 
obviously desirable that more accurate determinations should be made. 
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SUMMARY 

An investigation has been made of the friction and surface damage 
produced when a hard steel ball rolls between flat parallel surfaces of a 
softer metal. The behaviour may be divided into two parts. When rolling 
first commences, even at small loads, there is marked plastic flow and the 
resistance to rolling is due primarily to the plastic displacement of metal 
ahead of the ball. With successive traversals of the same track there is 
an increase in track width and a diminution in rolling resistance. An 
equilibrium state is ultimately reached where no further increase in track 
width occurs so that the deformations are essentially elastic. The rolling 
resistance is now constant and is unaffected by boundary lubricants. 
From a detailed study of the rolling of spheres and cylinders on rubber 
it is shown that there is no interfacial slip, the whole of the rolling resis- 
tance being due to elastic hysteresis losses. It is suggested that in the 
rolling experiments on metals within the elastic range the resistance to 
rolling is also due primarily to elastic hysteresis losses. Thus in rolling, 
there is no ‘ friction ’ in the conventional sense of the word. 


WE have recently carried out an investigation of the friction and surface 
deformation produced when a hard steel sphere rolls between flat parallel 
surfaces of a softer metal. The greater part of the earlier work was 
carried out by Mr. K. R. Eldredge on apparatus which he has already 
briefly described (Eldredge 1950). The behaviour may be divided into 
two parts. When rolling first commences, even at small loads, there is 
marked plastic flow and the resistance to rolling is due primarily to the 
plastic displacement of the metal ahead of the ball. Indeed, there is a 
fairly close quantitative agreement between the observed resistance and 
the product of the cross-sectional area of the grooved track and the 
yield pressure of the metal. In the first roll, the ball is supported on the 
front half of the circle of contact. On rolling back along this track, there 
is no such support for the ball, unless it sinks further into the metal. 
This leads to an increase in the width of the groove and further plastic 
displacement of metal. With successive rolls there is a slow increase in 
groove-width and a gradual decrease in rolling resistance. For a metal 
which does not work-harden appreciably (so that the area supporting 
the load is substantially constant for each roll) a simple theory may be 


ers Spee et ee eee 
* Communicated by the Author. 
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developed relating groove-width and friction as a function of the number 
of traverses. The results for tin are given in fig. 1 and it is seen that the 
agreement is reasonably good. Since this phase’ of the rolling process 
involves bulk deformation of the metal, the rolling resistance, as we 
might expect, depends little on the lubrication. With a metal like 
indium, however, the adhesion which is very marked for clean surfaces 
contributes some part of the rolling friction and this is naturally affecter 
resence of lubricants. 
ae Fallin proceeds, there is an increase in track-width as describec 
above, a slight change in track curvature and, accompanying thes 


Fig. | 
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Variation of track width and rolling resistance for repeated traversals of same 
track, Steel ball on tin surface. Vertical ordinates are arbitrary. 


changes, work-hardening of material around the track. As a result of 
these three factors a stage is gradually reached (after about 100 traverses) 
at which the whole of the load is borne by the ellipse of contact formed 
by elastic deformation within the width of the existing track. No 
further appreciable increase in track-width occurs and the rolling friction 
remains substantially constant for repeated traversals of the track. 
Here again, lubricants have practically no effect on the friction. 

In order to study rolling friction in the elastic range a more detailed 
investigation was made of the rolling of hard spheres and cylinders on 
rubber. More than 70 years ago, Osborne Reynolds (1875) pointed out 
that when a hard metal cylinder rolls on rubber it traverses a distance 
Jess than its periphery in a single revolution. This is because it measures 
out its periphery on the stretched portion of the rubber within the region 
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of contact, and Reynolds regarded this as equivalent to a small amount 
of slip between the roller and the rubber. He also provided an ingenious 
explanation of the independence of the rolling friction on lubrication. 
Recently, Dr. Ailsa Lowen, working in this laboratory, has found that 
although there is, as Reynolds observed, a difference between the 
circumference and the distance travelled forward in one revolution, there 
#s no slip between the rubber and the roller. The rubber is stretched 
vefore it enters and after it leaves the region of contact, and no inter- 
acial slip occurs. Similar effects are observed when hard spherical 
surfaces roll on rubber. Tangential forces are developed during rolling 
at the interface, but these are always appreciably less than the force 
“required to produce interfacial slip even in the presence of lubricant. 
This absence of slip between the rolling surfaces explains the fundamental 
observation that lubrication has no effect on rolling friction in the elastic 
range. (There is evidence that some adhesion may occur at the interface, 
so that some work may be expended in separating the ball from the 
surface as rolling proceeds. This will be reduced by lubricants but its 
contribution to the overall rolling resistance appears, in general, to be 
- negligible.) 

We must therefore seek a mechanism of rolling friction that does not 
involve surface slip or plastic flow. It is suggested here that the friction 
in the elastic range is primarily due to elastic hysteresis losses. If a 
steel ball, radius 7, rolls in a groove of radius r, under a load W, it forms 
an ellipse of contact of semi-major axis @ and semi-minor axis b, where 
the ratio a,/b) depends only on the ratio r,/r;. Using Hertz’s classical 
equations for elastic deformation (see Allan 1945) we have shown that the 
elastic work done in rolling forward a distance of one centimetre is 


On Weel 

p— T6 (=) [(#)- KH A/0o, A aed. eee) i (1) 
where f(£H)=(1—o,*)/H,+(1—o,”)/EH2, 01, o, being Poisson’s ratio : 
E,, E, Young’s modulus for the ball and groove material respectively, 
and & is an elliptic integral depending only on the ratio a,/by itself, 
that is, on the ratio r,/r,. For a ball rolling on a flat surface, the region 
of contact is a circle, so that ay=by) and 4 =7/2. 

If the hysteresis loss is constant, as it is for rubber at fairly constant 
rates of strain, the rolling friction F should be a constant fraction of @. 
Since a, is proportional to W1? it follows that Pa W**. Typical results 
plotted on a log-log scale for a steel ball rolling on flat rubber specimens 
of the same Young’s modulus but of widely differing hysteresis losses 
are shown in fig. 2. It is seen that F « W**, which is in good agreement 
with the theoretical slope. Further, direct determination of the 
hysteresis losses under normal loading, gave values of 11% for one rubber 
and 1-4% for the other. The observed frictional values agree with the 
calculated elastic energy values if hysteresis loss factors of 12% and 
1:3% respectively are assumed. The agreement is very satisfactory. 
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‘These ideas have been extended to metals when the track has reached 
its equilibrium state. For metals such as copper and nickel the ee 
friction is relatively high and falls off with load more rapidly than W248 ; 
indeed F is more nearly proportional to W?. Hysteresis losses can well 
explain these results if it is assumed (a) that for metals such as copper 
and nickel hysteresis losses of the order of 10% or more may occur, (b) 
that the hysteresis loss factor is not constant—as with rubber—but 
diminishes with decreasing strain. Both these assumptions are borne 
out by direct hysteresis experiments, on a standard Chevenard equipment, 
for metals deformed just beyond their elastic limit and then subjected 
to cyclic strains. Large hysteresis loops are observed which remain 
fairly constant in size for hundreds or thousands of cycles but fall off 


Fig. 2 
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Rolling resistance of 5/8in. steel ball rolling on rubber. x High hysteresis 
loss rubber. O Low hysteresis loss rubber. Note that the points are 
plotted on a log-log scale. : 


rapidly when the strain is reduced (Chevenard 1943). In this test the 
strains are torsional but the magnitude of the effect is confirmed. For 
typical ball-race steel the frictional resistance corresponds to a hysteresis 
loss of about 1%. This again is supported by direct hysteresis measure- 
ments on similar steels (Brophy 1936). If this value is used to calculate 
the frictional resistance of real ball or roller bearings it yields values of 
the right order of magnitude (ju 10-3). 

Although the behaviour of real ball and roller bearings is complicated 
by many other factors, the elastic hysteresis loss described in this com- 
munication must play a fundamental part. It is well known that ball 
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and roller bearings fail by sub-surface fatigue (Jones 1946) rather than by 
surface abrasion but the source of the rolling resistance itself has long 
been a matter of speculation. The work described here shows that in 
rolling there is no ‘friction’ in the conventional sense of the word. 
The resistance to rolling is primarily due to hysteresis losses. A fuller 
account of this work and a more detailed discussion of its practical 
implications are being published elsewhere. 
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ABSTRACT 


An anthracene scintillation counter has been used as a spectrometer 
to measure the energy distribution of initially monokinetic neutrons 
scattered from iron and aluminium. In each case discrete groups of 
neutrons with reduced energy were found, showing that the scattering 
nucleus was being raised to its excited states. The differential cross 
section for scattering through 97° has been measured for the various 


processes. 


§1. INTRODUCTION 


THE inelastic scattering of neutrons has been reported by many observers 
(Allen and Hurst 1940, Barschall and Ladenburg 1942, Barschall, Manley 
and Weisskopf 1947, Dunlap and Little 1941, Grahame and Seaborg 
1938, Mandeville and Swann 1951, Stelson and Preston 1952, Whitmore 
and Dennis 1951 and others) but only in comparatively few cases has. 
the method been capable of measuring the spectrum of the scattered 
neutrons (Dunlap 1941, Mandeville 1951, Stelson 1952 and Whitmore 
1951) and in only one case has the resolution been so good that individual 
groups were reported (Stelson 1952). 


§2. Tot NEvuTRON SPECTROMETER 


The chief difficulty in all inelastic scattering measurements is to find a 
satisfactory neutron spectrometer. During the following work an 
anthracene scintillation counter has been used in the manner already 
described by the author (Poole 1952). A cylindrical piece of anthracene 
lcm diameter by 0:5cm high was cemented to the window of an. 
E.M.I. type VX 5045 photomultiplier, and the pulse height spectrum 
arising from recoil protons in the anthracene was recorded with a 
30-channel pulse analyser. To convert this into the energy spectrum of 
the recoil protons the theoretical expression derived by Birks (1951) 
for the relation between pulse height and energy in anthracene was used. 
The result of a run allowing the 2-5 Mev neutrons from the D-D reaction 


* Communicated by the Author. 
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to fall directly on to the counter is shown in fig. 1. Figure 1 (a) shows 
the pulse height spectrum and fig. 1(b) the energy spectrum of recoil 
protons. The neutron spectrum is given in terms of the proton spectrum 
by the relation 


E d 


N(E)« G(H) de 


P(E), 


where N(£) is the neutron spectrum, P(Z#) the proton spectrum and o(Z) 
the n—p scattering cross section. The resolution at 2-5 Mev was found 
to be about +150 kev. 


Fig. 1 
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§3. SCATTERING GEOMETRY 


Before a scintillation counter of the type described can be used for 
scattering measurements, some method must be found to overcome its 
high sensitivity to y-rays. Hach neutron inelastically scattered leaves 
the target nucleus in an excited state which will usually decay to the 
ground state by y-emission. Thus the number of y-quanta from the 
scatterer itself will be of the same order of magnitude as the number of 
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scattered neutrons, and the counter sensitivity to, for example, a 1 Mev 
y-ray is approximately equal to its sensitivity to a 2-5 Mev neutron. 
In addition there will be a background of capture and scattering y-rays 
from all parts of the target room. 

The way the y-ray effect is eliminated is best seen from fig. 2. The 
target of the high-voltage machine is at A, and the D—D reaction was used 
to provide a monokinetic neutron source. As the target was thick, 
the scatterer (B) had to be placed in such a position that only the 
monokinetic neutrons emitted at 90° to the direction of the incident 
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deuteron beam were utilized, and it was therefore in the form of a ring 
placed round the target tube. Scattered neutrons were detected by the 
counter placed at K, directly below the neutron source, a polythene 
cylinder (D) preventing neutrons from reaching the counter directly 
from A. In addition either a 1-cm lead absorber (C, fig. 2 (a)) or an 
11 cm polythene absorber (C’, fig. 2 (b)) could be placed in the path of the 
neutrons scattered from B into E. A second scintillation counter was 
placed at F to give a count proportional to the neutron emission from A 
during any run. 

To obtain the distribution of recoil protons caused by scattered 
neutrons from B two runs were taken, one with the lead (C) in place and 
one with the polythene (C’). The difference between these is due only to 
neutrons from B if three conditions are satisfied. Firstly, the 
transparency of C and C’ to y-rays must be equal at all y-ray energies ; 
secondly, the general background of pulses caused by y-rays and neutrons 


Fig. 3 
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from the target room shall be unaltered when the lead and polythene are 
interchanged, and thirdly, the total neutron cross section for lead must 
be sensibly constant over the range of neutron energy concerned. 

The first condition cannot be satisfied exactly for all y-ray energies 
below 2-5 mev, largely because of the effect of photoelectric absorption 
in the lead. In order to assess the error from this cause, the ‘ effective’ 
sensitivity of the counter to y-rays when used in the way described has 
been calculated as a fraction of its direct. sensitivity. The y-ray 
absorption coefficients were taken from the work of Cowan (1948), and 
the lengths of the absorbers were chosen so as to make their transparencies 
exactly equal at 1:2mev. The result of this calculation showed the 
residual y-sensitivity to be 3% or less of the direct sensitivity in the range 
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0-8 mev to 3-0 Mev, rising to 5% at 0-6 mev. To check that interchanging 
the lead and paraffin absorbers did not alter the shape of the pulse-height 
distribution from y-rays a trial run was made with a ring of 6°Co sources 
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substituted for the scatterer B. This same run was used to adjust the 
length of the polythene absorber to give equal transparencies, and the 
result, after this adjustment had been made is shown Inet eres 
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The validity of the second assumption was checked by making runs 
without ring B in place, the results of which are plotted in fig. 4. The 
agreement of the two sets of points is satisfactory provided that no use 
is made of any points below 10 v bias (0-8 Mev energy). Finally, the 
third assumption was checked by making use of the figures published for 
the total cross section of lead (Adair 1950). This cross section varies 
smoothly from 5-4 barns at 0-8 mev neutron energy to 5-2 barns at 
2-5 Mev, giving a variation of neutron attenuation from 17-8% to 17-2% 
in the lead disc. 


§4. EXPERIMENTAL PROCEDURE 
For each run the following operations were carried out :— 

1. Counter E was placed at a known distance from the target, in a 
direction at 90° to the beam, and a pulse-analyser run made during the 
arrival of 5.10? monitor counts. This run served as an energy calibration 
for counter E and also related the sensitivity of the monitor (F) to that 
of E, which was needed to calculate the scattering cross section from 
the results. 

2. Counter, scatterer, polythene rod and lead absorber were set up 
as in fig. 2 (a) and a pulse analyser run taken during the arrival of 
5.104 monitor counts. 

3. The polythene absorber was substituted for lead, as in fig. 2 (b) and 
the run repeated. 

4, The channel widths of the pulse analyser were measured with a 
precision pulse generator. 

The last three operations were repeated about six times to reach the 
statistical accuracy required, and then the spectrum of recoil protons in 
the counter crystal, due to scattered neutrons from B (fig. 2) was obtained 
from the difference between 2 and 3 by correcting for the non-linear 
response of anthracene to protons. To derive the neutron spectrum two 
methods were used. In the case of runs made with an iron scatterer (fig. 5) 
a smooth curve was drawn through the recoil proton energy distribution, 
this curve differentiated and the neutron spectrum calculated from the 
differential coefficient. In the case of the aluminium runs (fig. 6) where 
the shape of the proton recoil curve was not so clear, it was felt that this 
procedure might lead to too much personal error, and an arithmetical 
method was used. Groups of five points were taken (~200 kev energy 
interval) and the gradient of the best straight line through these points 
determined by a least squares method. These values were then used to 
determine the neutron spectrum plotted in fig. 6 (6). This procedure 
corresponds to taking the mean gradient of protons over an. interval of 
+100 kev, each side of the point required, which has some justification 
as the resolution of the counter is known to be not better than +150 kev, 
and any fluctuations in gradient within this range are therefore 
legitimately ascribed to statistical deviations, and may be smoothed out. 


1066 M. J. Poole on the 


§5. RESULTS 
1. Iron Scatterer 
Figure 5 shows the proton recoil spectrum and the spectrum of scattered 


neutrons derived from this. Two groups of neutrons are clearly resolved, 
one at 2-35Mmev and the second at 1-55mev. The higher group 


Fig. 5 
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corresponds to elastic scattering of the incident neutrons (mean energy 
2-4 Mev), while that at 1-55 Mev may be attributed to inelastic scattering, 
leaving the iron nucleus excited to a level at 0:8+0-1 Mev. This is in 
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agreement with the work of Stelson and Preston (1952) who report a 
level at 0-85-+0-05 Mev, and with the measurements of Elliott and 
Deutsch (1943) who find a level at 0-845 Mev in Fe, by observations on 
the radioactivity of **Mn, and at 0-82 Mev by observations on Co. 


Fig. 6 
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2. Aluminium Scatterer 


Figure 6 shows the proton and neutron spectra for scattering from 
aluminium. There is a neutron group at 1-5 Mev, and a second group, 
ascribed to elastic scattering, at 2-3. Mev. Between these two is possibly 
a third group at 1-95 Mev. 
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The group at 1:5 Mev corresponds to an excitation of the residual 
nucleus of 0-9 Mev. This must be compared with the results of Reilley 
et al. (1952) on inelastic scattering of protons, giving levels at 0-81 and 
0-99 mev excitation, and also with the y-spectra of Hedgran, Hole and 
Benes (1948) and of Itoh (1941), who place the levels of 77Al at 0-84 Mevs 
and 1-85 ev. The level found here at 0-9 Mev is not in disagreement 
with this, particularly when it is remembered that two levels at 0-81 and 
0:99 mev would certainly not be resolved. The neutron group at 1-95 Mev, 
if present, would correspond to a level at 0-35 Mev excitation, which has 
not appeared yet in either the proton scattering or the y-spectrum work, 
but it is not sufficiently clearly resolved here to be certain that it is 
present. 


3. Cross Section Values 

The differential cross section for scattering 2-5 Mev neutrons through 
97° has been calculated for the various processes. These are tabulated 
below. 


Element Elastic Inelastic Comments 

Tron 0-10+1-01 barns/ 0-08--0-02 barns/ With 0-8 Mev 
steradian steradian energy loss 

Aluminium 0-13-+-0-03 barns/ 0-05--0-02 barns/ With 0-9 Mev 
steradian steradian _ energy loss. 


As the total cross section for iron at 2-5 Mev is 2-8 barns, and that for 
aluminium 2-5 barns, these figures seem high when it is remembered that 
a large proportion of scattering is commonly said to be potential scattering 
into the forward angles, but in the absence of any information concerning 
angular distributions, no definite conclusions can be drawn. 
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ABSTRACT 


The radioactive decay of 87Rb has been studied using the scintillations. 
from a crystal of rubidium iodide activated with thallium, some fluore- 
scent properties of which are described. This method free from source ~ 
absorption and scattering troubles, indicates a B-spectrum with an end 
point of 275 Kev. The Fermi plot gave a forbidden shape, for the most. 
part close to that given by recent proportional counter work using 27 
geometry, but having slightly more counts at the high energy end, and 
fewer at the low energy end; the half-life has been determined at. 
(5-90-0-3) x 1019 years. The theoretical interpretations are indicated. 


§1. INTRODUCTION 


THE exact nature of the decay process of **Rb has been in dispute. 
Conversion peaks were reported by Ollano (1941) and f-e coincidences by 
Haxel, Houtermans, and Kemmerich (1948). Curran, Dixon and Wilson 
(1951) using a proportional counter to determine f-ray shape found no 
evidence of conversion lines, and McGregor and Wiedenbeck (1952) have 
recently shown electron—electron coincidences unlikely. The end point 
which was thought by Libby and Lee (1939) to be 130 Kev has been 
more recently specified at 270-275 Kev (cf. Bell, Cassidy and Davis 
1950, and Curran, Dixon and Wilson 1951). Several estimates of life- 
time have been made, generally dependent on correction factors ; and 
recent values lie in the neighbourhood of 6 x 101° years. 

In view of the various results, and considering the methods to date, 
it was felt that an investigation using an activated crystal of rubidium 
iodide would help in several ways, providing suitable scintillations could 
be obtained. The absence of self-absorption and back-scattering errors 
would make this 47-method especially suitable for estimating lower 
energy particles without distortion, and also for yielding an estimate of 
lifetime. Further the higher energy particles would be counted fairly 


completely in a crystal; and the radiation background could be made 
negligible. 


§2. GENERAL Mernop, PROPERTIES OF THE CRYSTALS AND CALIBRATION 


A few grams of rubidium iodide activated with 0-5% thallium iodide 
were placed in a quartz container which was evacuated near the melting 
point of the mixture and sealed. The whole was made to traverse a 


ee ee 
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temperature gradient which permitted melting and slow cooling to room 
temperature. ‘Transparent, glassy crystals were obtained, relatively free 
from hygroscopic tendencies. 

Such a crystal was mounted in paraffin on an E.M.I. multiplier type 
VX 5032 and covered by an aluminium reflector. The system was 
enclosed by a two-inch thickness of lead. The amplified pulses led into 
a single channel ‘kicksorter’, activated by a gating unit normally 
triggered by the pulse itself. The pulses observed under y-ray bombard- 
ment were of height ¢ those of thallium activated sodium iodide, and had 
a rise time better than about 1 sec. It was early observed however 
that a phosphorescence decay was also produced lasting the order of a 
millisecond. Crystals prepared with considerably less than 0-5°%, and 
with 2% of thallium iodide also showed this phosphorescence. To 
investigate and avoid this radiation, a paralysis circuit, fed by a discri- 
minator, could be used to operate the gating unit of the kicksorter. 
It was found that the phosphorescence had no ill effect except in studying 
the very low energies, where it gave, in the absence of paralysis, spurious 
counts. At the lowest energies paralysis of order 1-3 milliseconds enabled 
the investigation to be carried out to regions where the tube noise became 
a final barrier. 

The main pulse was shown to vary linearly with energy for electrons 
moving in the crystal, over a range 17 Kev—630 Kev, using various X- 
and y-rays and higher energy electrons. Calibration was effected using the 
fluorescent Kax-rays of molybdenum, the fluorescent Kx-rays of silver, 
114Tn (Kx and conversion electrons), RaD (46-7 Kev y-ray), "Hg (Kx and 
y-ray), 1°*Cs (conversion electrons). Calibration with x- and y-rays 
generally involves production of a photoelectron and subsequent capture 
of the associated x-ray. This would imply linearity of calibration here 
some Kev below the figure of 17 mentioned above. 


§3. REsuLTS 

The pulse height distribution using a crystal 4mm 4mm x2} mm 
is shown in fig. 1. The curve is composed of three sets of points. The 
amplifier gain was greater for the two sets at lower energy. ‘The points 
indicated by triangles were obtained using paralysis—other points (omitted 
for clarity of drawing) at higher energies using paralysis fell on the curve 
drawn. Reliable calibration was not readily possible below ~17 Kev, 
but one point has been included in fig. 1, below this (viz., at ~13 Kev) 
assuming the linear energy relation to continue. 

The curve of fig. 1 reaches background just after 300 Kev, but the 
character of the curve appears to alter at about 250 Kev, indicating an 
end point near 275 kev. The straggle here is due to the resolving power 
of the crystal—the *°Hg y-ray peak at 280 Kev had half width at half 
height ~15%. The value of 275 Kev is confirmed from the Fermi 
plot below. The curve is smooth in form and further intermediate 
points have been taken which support this view. This confirms the 


1072 G. M. Lewis on the 


absence of conversion peaks and y-rays in the range investigated. With 
47 geometry, these would appear as part of an integrated spectrum 
(cf. Bannerman, Lewis and Curran 1951). Preliminary work with another 
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crystal gave similar results. It should be pointed out here that Curran, 
Dixon and Wilson (1951) could find no evidence of K and L X-rays in 
their investigation. 
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The energy distribution N(W) dW can be written :— 
N(W) dW=CF(W, Z)(Wo—W)(W2—1)!?W dW—(1), 
where Wo is the maximum energy and F(W, Z) is the Coulomb correction 


factor; C is constant for allowed transitions, but is more generally, 
energy dependent. The Fermi plot given by 


(V(W)/F(W, Z)(W?—1)3? W 3? 
is shown in fig. 2, where the Bethe-Bacher approximation has been used 
for F(W, Z). Points in the straggle area have been omitted. The Fermi 
plot indicates an end point of 275 Kev. The curve shows good agreement 
over most of the energy range with the proportional counter curve shown 


Fig. 2 
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Solid curve —Fermi plot of ®’Rb B-spectrum obtained with activated rubidium 
iodide crystal. .—.—. Fermi Plot obtained by Curran, Dixon and 
Wilson (1951) with the proportional counter. 


for comparison in fig. 2. There are slightly more counts at the higher 
energy. This is to be expected as these particles are absorbed fairly 
completely in the crystal and not so well in the gas counter. There are 
less counts at the lower energy. This is due to the absence of back 
scattering, and also to the lower resolution of the crystal method at these 
energies. 
The total count above 17-4 Kev has been determined using various 
paralysis times in the neighbourhood of 0-3 milliseconds. Allowing for 


1074 On the Natural Radioactivity of Rubidium 


loss of proper counts through coincident arrival within the paralysis 
time, the crystals of weight 0-1487 gm gave a total count above this 
energy of 1975/min. To estimate the overall counts, extrapolation of 
the Fermi plot along the dotted line to zero energy 1s necessary, involving 
an additional 600 counts/min. A half life figure of (5-90+-0-3)10 years 
is obtained when the 87Rb abundance figure is taken as 27-2%. 


§4. DIscuSssION 


With an end point energy of 275 Kev and half life of 5-90 x 10" years, 
log ft.=17-5, and the transition **Rb— ‘’Sr is third forbidden. The 
nuclear spin change between the ground states, 3/2 9/2, accords with 
this. According to the nuclear shell model, the transition, p3/2 > g9/2, 
involves a parity change. Theoretical correction factors (involving C' of 
eqn. (1)) to the Fermi plot shape have been formulated by Greuling 
(1942). Pseudoscalar interaction (3P) is eliminated on parity considera- 
tions. Calculation shows that scalar interaction (3S) does not explain 
the 8’Rb plot either. Generally the correction factors involve knowledge 
of certain nuclear matrix elements ; Tomozawa, Umezawa and Nakamura 
(1952) have recently shown that this type of plot for *’Rb could be given 
by and only by polar vector (3V) or tensor (37') interaction. Their 
theoretical data, involving one chosen parameter, has been derived for the 
region 50-280 Kev, and gives good fits with both the experimental curves 
of fig. 2; thus division of the ordinates of points on these curves by the 
root of their correction factor (1/C) gives approximately linear plots. 
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SUMMARY 


The deformation of silver wires at high temperatures and low applied 
stresses has been investigated experimentally. Considerable evidence has 
been obtained that the crystal grains in the wires deform by the self- 
diffusion mechanism postulated by Nabarro and considered in detail by 
Herring (1950). A new explanation of the development of kinking is 
advanced to account for the experimental observations. It is shown that 
surface energy effects can give rise to offsetting in a direction which causes 
the wires to shorten. 


$1. IvtRopUCcTION 


In a recent theoretical paper on the deformation of polycrystalline solids 
at: high temperatures under small stresses, Herring (1950) considered in 
detail the case of wires in which the crystal grains occupied the entire 
cross-section. He suggested that in this case the grains themselves would 
deform by the self-diffusion mechanism proposed by Nabarro (1948), the 
migration of vacant sites through the crystal lattice effectively transporting 
material in the opposite direction. Assuming that the external surfaces 
and grain boundaries could act as sources or sinks of vacant lattice sites, 
he derived an expression for the magnitude of this effect which compared 
favourably with existing experimental observations. 

In addition, he predicted that offsetting and kinking of the wires would 
take place under certain conditions. Offsetting is produced when slip 
occurs in the plane of a grain boundary whose normal is inclined to the 
wire axis. Herring considered that this effect would be observed only 
when the applied stress was somewhat larger than that necessary to prevent 
the shortening of a long grain under its own surface energy. inking, or 
the development of an angle between the axes of adjacent grains in the 
wires, could be produced in two ways. A relative rotation of two grains 
about an axis normal to the wire axis would give rise to kinking, but 


_ Herring considered that this process would be slow. Secondly, the 


anisotropy of grain boundary energy could cause the relative rotation of 
two grains about a normal to a plane grain boundary inclined to the wire 
axis. This would produce a kink, and also an axial twisting of the wire, 
but since no transport of matter would be required, the rate at which the 
kink developed might be faster than that of the first process. 
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The deformation of copper wires at high temperatures has been studied 
by Udin, Shaler and Wulff (1949). Alexander, Dawson and Kling (1951), 
working with gold, reported a,lack of reproducibility of strain measure- 
ments on their wires, which they attributed to kinking, to offsetting, and to 
the necking which developed in some of their wires. Recently, the 
surface ‘ tension ’ of the solid metals has been estimated from observations 
on wires by Funk, Udin and Wulff (1951) for silver, and by Buttner, 
Udin and Wulff (1951) for gold. Buttner, Funk and Udin (1952) investi- 
gated the mechanism of deformation of gold wires, and showed that the 
diffusion mechanism would account for the creep observed at stresses less 
than about 5-5 10° dynes/cm?, but that above this stress an additional 
slip mechanism occurs. They reported the development of kinks and 
offsets in their wires, but failed to observe the bulging and necking 
reported by Alexander, Dawson and Kling, and observed no axial rotation 
of the wires at the grain boundaries. 


§2, EXPERIMENTS ON POLYCRYSTALLINE SPECIMENS 


In these experiments, the specimens resembled the copper specimens of 
Udin, Shaler and Wulff. They were prepared from cold-drawn silver 
wire of 42s.w.g. (0-01016 cm diam.) supplied by Johnson Matthey and — 
Co. Ltd., who estimated the purity to be greater than 99-999%. A knot 
was tied in the wire at one end, and the wire cut at a measured distance 
from this knot. A second knot was tied about 10 cm below the first, and 
the remainder of the wire coiled to form a weight. The knots were 
arranged to be as nearly as possible in the same plane, so that any appreci- 
able twisting of the wire about its axis could be observed. 

After being straightened as far as possible by hand the specimens were 
hung on silver supports and lowered into a vertical silica tube furnace in 
which an atmosphere of oxygen-free nitrogen, purified as already described 
(Greenough and King 1951), was then established. The specimens were 
first annealed for 1} hours at the temperature of test. When cool, the 
wires were raised out of the furnace and examined, the gauge length being 
measured by a travelling microscope, and the relative planes of the knots 


noted. These observations were repeated at the end of each successive 
annealing treatment. 


(a) Kinking 


At the end of the preliminary annealing treatment the wires were quite 
straight. Subsequently kinks developed, and became apparent in all 
specimens when the strain exceeded about +0-1%, but no appreciable 
axial twisting of the specimens was observed during this process. 

It was found that the grain boundaries in the wire could be seen 
because of the thermal etching which took place in the furnace. If the 
wire was mounted on a glass slide and examined microscopically with a 
4mm objective and the illumination coming from behind, enough light 
was reflected by the surface of the objective to enable the grain structure 
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of the wire to be seen. Simultaneously the outline of the wire was seen 
against bright ground (figs. 2, 3, 4 and 7, Plates LXIX and LXX). 
Although the grains generally occupied the whole cross-section of the 
Wire, a small number occupied only a part of the cross-section, producing a 
forked grain boundary as shown in figs. 2and 3. In these cases, there are 
two grain boundaries on one side of the specimen, but only one on the 
other. Thus one side of the specimen will tend to deform at twice the 
rate of the other if the grain boundaries are acting as sources or sinks of 
the diffusion currents, and a kink will tend to develop. In wires which are 
shortening under surface energy effects, the small grain will be on the 
concave side of the kink (fig. 2), but on the convex side if the wire is 
increasing in length (fig. 3). The stresses set up in the wire below the kink 
could cause kinking in the opposite direction to take place at plane grain 
boundaries as shown at Bin fig. 1. In all cases, the angle of kink is limited. 


‘Rigged Fig. 6 
<A 
<A 
<—B 
x 50 . x 50 
After 85 hours at 920° oc. After 61 hours at 920° c. 
Wire shortening. Wire extending. 


~ When the wire is contracting, the equilibrium angle of kink is reached when 
the reduction in energy due to the destruction of the wire surface on the 
inside of the kink is balanced by the increase in energy due to the raising 
of the centre of gravity of the wire below. When the wire is extending 
under the applied load, energy is required to create the extra surface on the 
outside of the kink and this must come from the lowering of the centre of 
gravity of the wire below, which occurs only when the angle of kink is 
small. It appears that the kinks in gold wires illustrated by Alexander, 
Dawson and Kling can be explained on the basis of this theory. As no 
axial twisting of the wires was noticed during the development of the 
kinks, it seems unlikely that either of the rotation mechanisms postulated 
by Herring leads to the development of kinks to any appreciable extent. 
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(b) Offsetting 

The offset shown in fig. 4 cannot have developed by the mechanism 
postulated by Herring since the direction of offsetting is such as to re 
the weight of material under the plane of slip. However, slip must ave 
taken place since the centre-line of the wire is offset. Consider a grain 
boundary inclined to the axis of the wire, shown in section in fig. 5 (a). 
At equilibrium in the absence of an applied stress, the plane of the grain 
boundary bisects the angle between the two external surfaces (Chalmers, 
King and Shuttleworth 1948, Greenough and King 1951). Thus equili- 
brium could be attained either by the development of bulges (fig. 5 (b)), or 
by slip taking place in the plane of the boundary, as indicated in fig. 5 (c). 
Of these two processes, the latter results in a final state of lower energy than 
the former, since the resultant increase in external surface area 1s Jess and 
the grain boundary area has decreased. This is the situation in fig. 4 
where the load on the wire was very small compared with the surface 


energy effects. 


Fig. 5 
(a) (0) (c) 
No equilibrium. Equilibrium. No offsetting. Equilibrium. Offsetting. 


As the tensile stress in the wire increases, there will be an increasing 
tendency for slip on the plane of the boundary to increase the length of 
the wire (figs. 6 and 7), until, as Herring points out, the tensile stress is 
sufficient to provide the energy necessary for the grain boundary to become 
two external surfaces. At a stress of this magnitude, or greater, no 
equilibrium is possible and the two grains could slide apart at the grain 
boundary. 


(c) Necking 

Wires which broke under the applied stress usually showed a series of 
necks near the fracture. Each neck occurred at a grain boundary, as 
would be expected if the grain boundaries were serving as sources for the 
vacant lattice sites. A crystallographic slip mechanism (Cottrell and 


Jaswon 1949) would be expected to lead to necking at the centre of the 
grains. No bulging was observed in any of the specimens, 
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(d) Rate of Strain | 


For a wire of radius r, and plane grain boundaries normal to the axis 
of the wire spaced a distance L apart, Herring derived an expression (his 
eqn. 16) from which the rate of strain under a given stress may be 
calculated. Using the data of Hoffman and Turnbull (1901) 5 for 
the self-diffusion coefficient of silver, this reduces to : 


Strain rate=2-05 x 10-18 S/Lr, per sec, ... . . (1) 


for silver wire at 920°c, when L>2r, and S is the applied stress, corrected 
for surface energy effects (Udin 1951), in dynes/em?. 

Much quantitative data was obtained on the strains shown by wires, but 
the reproducibility of the results was poor. Specimens heated simul- 
taneously, and the same specimen on two successive heats, would show 
very different strain rates. The values of the surface energy of silver 
given by successive experiments with the same specimens steadily 
increased so that the final value would be two or three times the original. 
These effects may have been due to temperature variations, since an 
increase in temperature of 10°C increases the strain rate by about 17%. 
There was evidence to suggest that specimens heated simultaneously were 
not all at the same temperature, and it was known that the temperature 
of the furnace was not constant over long periods. However, it seems 
unlikely that this was the main cause, for with very good temperature 
control, the results obtained for gold by Buttner, Funk and Udin also 
show poor reproducibility, as can be seen from their figs. 2 and 3. 

It is nevertheless interesting to note that the average strain rates 
observed in the silver wires agreed with those predicted by Herring. 
Table 1 gives the strain measurements on nine similar specimens. Other 
experiments indicated that the surface energy of silver at 920°C was 
1 500 ergs/cm?+50%, a value in agreement with that obtained by Funk, 
Udin and Wulff. Thus, in eqn. (1), S is —22-0x 10* dynes/cm?, and the 
calculated strain rate is —3-00 x 10~® per sec, compared with the observed 
rate of —7-61x10~-® per sec. In a second experiment the calculated 
strain rate was —5:0210-® per sec, and the observed strain rate 
 —10-6x 10-° per sec. 

The agreement between the experimentally observed and the theoretic- 
ally predicted strain rates is satisfactory, when the accuracy of the experi- 
mental data and the nature of the theoretical calculations are cousidered. 
The expression used to calculate strain rate takes no account of kinking, 
offsetting, and the departure of some of the grain boundary shapes from 
the ideal transverse planes, all factors which would increase the calculated 

strain rate. 

The strain rate observed in copper wires by Udin, Shaler and Wulff 
was greater than the calculated strain rate by a factor of four to ten. In 
gold wires, the observed strain rate was about 30% less than the calculated 
strain rate (Buttner, Funk and Udin 1952). 
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§3. EXPERIMENT ON SINGLE CRYSTALS 


: An experiment was carried out with two single crystal wires of assay 
silver, 0-1 cm diameter, grown from the melt in graphite moulds by 
the travelling furnace method. For comparigon, a specimen which 
recrystallized to give several grains in the cross-section was made from 
hard-drawn assay silver wire of the same diameter. Silver weights 
hung from hooks formed at the lower ends of the specimens were used 
to provide loads. The gauge length was taken as the distance between 
two loops of fine wire sintered on to the specimens, this method being 
adopted to avoid recrystallization in the gauge length of the single 
crystal specimens. 

Table 2 compares the behaviour of the two single crystal wires with the 
polycrystalline specimen. The three specimens were annealed simul- 
taneously for 282-5 hours at 910° c. 


Table 2 


Applied stress at | Initial gauge} Increase of 
middle of gauge length gauge length 
length dynes/cm? cm cm 


Single crystal 1 131 000 4-044 0-005 
Single crystal 2 114 000 4-339 0-001 
Polycrystalline wire 126 000 4-492 + 0-058 


The maximum possible error in the measurement of the change in 
gauge length was +0:012cm. The applied stress will be reduced by 
approximately 30 000 dynes/cm? in each case by surface energy effects. 
The effect of the grain-boundary energy in the case of the polycrystalline 
wire will be small. Herring has pointed out (private communication) 
that the observed creep rate in the polycrystalline wire would agree with 
_ the theoretical rate (Herring 1950, eqn. 11) if the average grain diameter 

were of the order of a fifth of the diameter of the wire. This was 
experimentally confirmed. 

The stress applied to the single crystals was considerably less than that 
which would be expected to produce appreciable micro-creep (Chalmers 
1936, 1937). A stress of about 10’ dynes/cm® gave a strain rate of about 
5x 10-® per sec in tin single crystals. 


$4. CONCLUSIONS 


An experimental study of the high temperature deformation of silver 
wires under low stresses indicated that the principal mechanisms which 
operate are as follows : 

(1) The crystals themselves deform, probably by the self-diffusion 
mechanism first suggested by Nabarro. The diffusion currents flow 
between the external surfaces and the grain boundaries. ‘This was shown 
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experimentally by the behaviour of single crystals, the mode of develop- 
ment of kinks and necks, and the agreement between the experimentally 
observed strain rate and the strain rate calculated by Herring. It seems 
unlikely that dislocations in the crystal play any appreciable part in the 
diffusion process under the conditions of the experiment. 

(2) Deformation by offsetting takes place at grain boundaries inclined 
to the wire axis in such a way that the overall length of the wire either 
decreases or increases, depending on the relation between the applied 
load and the surface energy. 

(3) It is unlikely that kinking is produced by the mechanisms suggested 
by Herring. Kinking was observed, however, and an explanation of this 
is proposed. 
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SUMMARY 


The calculations of Montroll of the moments of the frequency spectrum 
of a number of crystal lattices enable an accurate estimate to be made 
of the zero point energy of these lattices. It has been customary to use 
the value 340, based on the Debye approximation, to estimate the zero 
point energy, and it is shown that this formula is only accurate if @,,, 
the limiting value of @ for high temperatures is used. It is shown that 
@,, can be simply expressed in terms of the force constants between 
atoms of the crystal. 


§1. THE Zero Point Enerey 
THERE are a number of substances in which zero point energy plays an 
important part in the behaviour of the solid state at low temperatures, 
and recent experiments by Dugdale and Simon (1952) have shown that 
for solid helium under pressure zero point energy plays a dominant part 
even at relatively high temperatures. In calculating the value of the 
zero point energy it has been customary to use the Debye approximation 
for the distribution of normal modes in the solid (Simon 1934, London 
1936). A normal mode of frequency v contributes $hv to the zero point 
energy, so that the total zero point energy per atom of a crystal is 
Le ldD a 

y8 iy hv? dv= gvp= 5k0. Memeo ees Pray: wb) 

The Debye approximation takes correct account of two important 
features of the true spectrum, the distribution for sufficiently low 


dy 


- frequencies, and the total number of normal frequencies of the crystal. 


However, as a detailed representation of the spectrum over the whole 
frequency range it is not very adequate (Blackman 1941), and, since low 
frequencies make almost no contribution to the zero point energy, it would 
be rather surprising if the formula (1) were at all accurate. 

In fact there is little difficulty in obtaining a much closer estimate of 
the value of the zero point energy if we make use of some results of 
Montroll (1942, 1943, 1944) on the nature of the true spectrum g(v) dv. 
Montroll employed a simple analytic method first introduced by Thirring 
(1913, 1914) for deriving the even moments of the spectrum, 


1 (ee) 
bax gy |, vale) de (even): = oem (2) 


* Communicated by the Authors. 
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and he applied this to the determination of the first seven moments for 
several different lattices. We shall use throughout the closely related 


moments, 

Hiy= fate) de= i (=, fee)=raalint))- « « @) 
0 Vm Vm : 
(Note that the parameter v,, appearing in these expressions is not the 
Debye vp, which is artificially defined in such a manner that integration 
of the Debye v? distribution over all v gives 3N ; on the more exact 
lattice theory v,, is the true maximum frequency, defined in terms of 
the force constants.) Essentially the zero point energy is given by the 
first odd moment, i.e. E,=3hv,,I(1), and the problem reduces to esti- 
mating J(1) from the known values of I(k). We have made use of two 

independent methods for this estimation. 
(a) The origin is transferred to v,,, and the modified moments, J(k), 
are calculated from 


Hk)=[ 2"*f(e) de (2781122), 


Then w can be expanded as a power series in x’?, 
¢=(l—2'2)12=1—fy'2—_ 44). 


The J(k) are all positive and decrease rapidly because of the behaviour 
of g(v) near v=0; we can thus substitute the known values of J(k) and 
obtain the first few terms of a rapidly converging series for J(1). The 
terms after the first are consistent in sign, and it is possible to estimate 
fairly accurately the contribution of the higher terms. This method is 
somewhat tedious, since the calculations must be performed afresh for 
each lattice, but it is susceptible of high accuracy. 

(b) (Suggested by Professor M. H. L. Pryce.) The value of J(1) is 
determined by a Lagrange interpolation from the known values of I(k). 
From the known behaviour of g(v) near v=0 it is clear that J(k) will have 
an infinity at k=—3; for the purpose of interpolation it is therefore 
better to use (k+-3) I(k). The calculation involved in this method is 
considerably simpler than in the previous method, and once the Lagrangian 
coefficients have been determined, they can be used for all lattices. 

Both methods were employed for the calculation of the zero point 
energy of the face centred cubic lattice, taking only nearest neighbour 

interactions into account. The results differed by less than 1 part in 
* 1000 (as will be seen from the next section), and it therefore seems that 
one can apply method (b) to other lattices with considerable confidence. 


§2. DrramEeD CaLCULATIONS FOR VARIOUS LaTTicEs 
‘ Montroll (1943, 1944) has derived moments for the simple cubic and 
body centred cubic lattices taking into account first and second neighbour 
interactions. If in a simple cubic lattice first neighbour interactions 
are zero, the lattice breaks up into a pair of independent face centred 
cubic lattices. We can thus from Montroll’s calculations deduce the 
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moments for the face centred cubic lattice with nearest neighbour 
interactions only, and they are as follows : 


Ho=1, fe=48, oy=208?, g=1148%, e—7038%, 
H19= 454985, 242=308538%, 44206757387, Pe cane) 
(B=«a/47?m). 
Here « is the force constant between nearest neighbours, m is the particle 
Mass, v,,—=1/74/(2«/m). 
Applying method (a) we derive the following series for the zero point 
energy : 
y= 3hr,{1—0-25—0-039063—0-013428 
— 0:006094— 0-003220—0-001880—0-001180...] . . . (7) 
The terms are decreasing steadily, and can be fitted asymptotically to an 
expression of the form A/n’. Extrapolating in this manner, we estimate 
the sum of the series to be 
Ey=1-0225hy,,. Merete re est (8) 
The terms involved in the extrapolation amount only to 0-0050, so that 
the correction is quite a small one. 
Method (b) gives J(1) as a linear combination of the known I(k) in 
the form 
I(1)=0-209473-+ 1-466309/ (2)—1-4663091(4) 
+ 1:4663091(6)—1-0473631(8)-+ 0-4887701(10) 
SO tess0ii (La) O-ONOLISI(I4)o6 1. Ge. (9) 
This formula is applicable to all the three dimensional lattices for which 
the moments have been determined, and in the particular case under 
consideration gives 
Pres (2002 0h wei Osa) ho fide. 10) 
It is interesting to check the interpolation method for the simple cubic 
lattice with nearest neighbour interactions only. The spectrum is in 
this case equivalent to that of a one dimensional lattice, so that the 
zero point energy can be calculated exactly. We should use (+ 1)I(‘) 
in the interpolation formula, and we derive the estimate 
Beebo ihyee ee. ere... (il) 
whereas the correct value should be 


3 
Ly= = hr, =0-95493hy,,, tive eee te ea) 


§3. Tue Usz or 0, 


The detailed behaviour of @ as a function of temperature for a face 
centred cubic lattice has been calculated by Leighton (1948). He finds 
that after an initial variation of about 15% in the range 7’<@/10, O 
settles down to a steady value; we shall denote by @,, the theoretical 
limiting value of @ as 7’ ©, and for the face centred cubic lattice 
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@ rapidly approaches’ @,, when 7’>6/5. The calculations of Blackman 
(1937) for the simple cubic lattice, and of Smith (1948) for the diamond 
lattice seems to indicate that this behaviour is fairly general. 

There are anumber of properties of @,, which are of considerable practical 
benefit. @,, can be defined even for systems whose spectrum at low 
frequencies does not conform to the Debye law (e.g. a simple cubic lattice 
with nearest neighbour interactions). It will be shown shortly that 
@,, can be very simply calculated for any lattice from given force constants, 
and that the formula (1) for the zero point energy provides a good 
approximation if @,, is used. 

At high temperatures the specific heat can be expanded as a series 
(Thirring 1913) : 


Cre Mofh\? | ty [h\4 Me (h\8 
sye—!— 3(a) +55 () — ots (a) meet 


h \2n : 
a) , the B,, being 


the Bernoulli numbers B,=%, B,=3, Bs=%, giving in terms of the I(k) 


where the general terms is given by (—1)" ral fen ( 


cs 7 te) hm\? , L(4) (hvm\* — £(6) (hrm\® ; 
svk—'— a2 (er) * 2a0\er) — eoas\er) ~ * (8) 
For the Debye spectrum the moments are given by 
3p" 
Me= TB Ip(k)=3/k+3. ee ee en 


Equations (15) and (14) can be used to derive an expansion of © in 
inverse powers of the temperature, and in particular, retaining only the 
first term in the limit as 7’—> 00, we find that 


hh fdpy  hvy [51 (2) 
ep k 5 = = at ei 3 . . . . . . (15) 


Formula (15) furnishes a simple means of calculating @,, for any lattice 
with given force constants. jr, can easily be deduced from the original 
panaeons of motion of the lattice. If these are written in the matrix 
orm 


Az=4r*my*z,: \.* ase event ee 


ly 1s essentially the diagonal sum or trace of A. It should be observed 
that the matrix A entering here is of order 3N, and is not the same as 
the matrices used by Thirring, Blackman and Montroll which are of order 
3 (or 38s when there are s atoms in unit cell). This matrix can be used 
to derive the first few moments without appeal to the cyclic boundary 
pera the complication increases very rapidly, however, for the 
Goquetey, is irlevat to the disbusiod: aia ates Rea ROE 
to facilitate computation and numeral! soeanny Tee 
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For the face centred cubic lattice considered in § 2, [(2)=4, and we 
find that ¢ 


and this is in close agreement with the result obtained by Leighton (1948). 
This can be written, substituting for v,,, or using j1y directly : 


ae wT °K 2) 


0-9129h /2a\1/2 
en) 

The numerical factor 0-9129, or \/%, occurring in eqn. (17), would seem 
to apply more generally since [(2) is found to be equal to 4 for each of 
the three cubic lattice structures, irrespective of ratio of first to second 
neighbour interaction. This means that v,, is related to the trace of the 
matrix for the lattice in each case by the same numerical factor; the 
reason for this simple relationship would appear to merit further 
investigation. 

We could equally well take account of distant neighbour interactions 


and write more generally 
h (Xa ;\12 
@,,=0-4109 7 (=) Aare te gy otoes - Pre (19) 


where «, is the ‘force constant ’ of the ith neighbour interaction which 
can be readily determined from the equations of motion. The lattice 
sums involved in the «,; have been evaluated for a number of lattices by 
Lennard-Jones and Ingham (1924) for potentials involving power laws. 
For an attractive law falling off as the inverse sixth power of the distance, 
the «, fall off as the eighth power, and the effect of distant neighbours 
on @,, in the face centred lattice is small (~3%). 

The above method of calculating @,, (and thus the zero point energy), 
is substantially simpler, and probably more accurate, than previous 
methods which made use of an average over elastic constants (Herzfeld 
and Goeppert-Mayer 1934, de Boer and Blaisse 1948). 

The method would also seem to be applicable to mixed crystals 
containing different atoms, using the matrix of eqn. (16). When these 
atoms form a well ordered mixture, the methods of lattice dynamics 
can be used, although the spectrum becomes more complicated (Smith 
1948). But when the mixtures are random, the symmetry of the system 
disappears, and it is difficult to see how to use classical theory without’ 
a large increase in labour. The second moment can still be determined 
very simply from the trace of the above matrix, and hence @,, can be 
calculated. If, for example, we assume that the force constants between 
all atoms are the same, but only their masses differ, we should have 
instead of (19), he 

6,,=0-4109 5 (3x,)"" ies ) pe ere (90) 


M4, Mo 


1088 C. Domb and L. Salter on the 


p, and p, being the fractions of the two constituents. It is hoped to 
apply this method to the effect of vibration terms in the theory of regular 
solutions, a problem which has been discussed recently by other writers 
(Prigogine and Garikian 1950, Rushbrooke and Ziman 1952). 


§4, Tue 340 FormuLA 


If we use the values of @,, given by (17), and estimate the zero point 
energy from the Debye formula (1) we obtain the value Ey=1-0270h ry, 
It will be seen that this is close to the value (8) given by refined calculations. 
If, instead, we had used the value of @ at 7’=0 given by Leighton’s 
calculations, we should have obtained a result in error by 15%. It is 
fairly clear that the variation in @ in the neighbourhood of 7’=0 is 
caused by the long wave region of the spectrum, and we should not 
expect this to contribute much to the zero point energy. 


Table 1 
Ratio Zero Point/hvy, 


Lattice 2nd : Ist Energy 2h0 n/hvm °% Difference 
Neighbour | (Interpolation) 


S.C. Py 1:0344 1-0270 0-7 
S.C. 1 1-0329 1-0270 0-6 
SHO 2 1:0297 1-0270 0-3 
B.C.C. 0 0:9956 1-0270 3-1 
B.C.C. i 1:0242 1:0270 0:3 
B.C.C. 1 1:0276 1:0270 0-1 
LEXORG. 2 1-:0226 1-0270 0-4. 
F.C.C 0 1:0220 1:0270 0:5 
S.C. 
(Linear Chain) - 0 0-9549 1-0270 ree) 
[3] Einstein 
Oscillator 0 4524 3:2 


1-5000 1-45: 


An alternative argument favouring the use of @,, for estimating zero 
point energy is as follows: by our determination of @,, we have ensured 
that the moments (0) and J(2) for the correct distribution and equivalent 
Debye distribution are the same, and we should not therefore expect (1) 
to differ very much. As a test of the formula we have compared the 
results with those derived from the interpolation formula (9) for a number 
of particular cases in which the moments had already been derived 
(Montroll 1943, 1944). The results are collected in table 1. The final 
two entries correspond to spectra which do not follow the Debye law, but 
even here it will be seen that the error is surprisingly small. In the other 
cubic lattice cases the only error greater than 1°% corresponds to the 
body centred lattice with no second neighbour interaction. For this 
lattice, however, the second neighbours are not much further distant 


than the nearest neighbours, and their interactions should be taken into 
account in practice. 
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§5. High TEMPERATURE VARIATION OF O 


The small variation in © at high temperatures was referred to at the 
beginning of §3. Formula (13) can be used to derive an expansion of 
@ in inverse powers of the temperature, and the first two terms are 


) 1 5 257 | (hen\? 


The expression for the higher terms rapidly becomes very complicated, 
but in the particular case of the face centred cubic lattice we have - 
evaluated the first seven terms, as follows: 


@,.=0-9129x[1—7-441 x 10-4(y/T)2-+ 1-447 x 10-5(y/T’)4 
— 1-549 X 10-7(x/T')®§+ 1-735 x 10-%x/T)8 
— 1-463 x 10-11(y/7) 94. 4-928 x 10-13(y/T)2 ... J. (y=hv,,/h) 
(22) 


It will be noted that the coefficients of the inverse powers of 7’ are 
extremely small. The formula can be used to evaluate @ from T= c 
to 7’ % 9/6 and it is found that in the whole of this range the variation 
is less than 2%. These results are in agreement with the calculations 
of Leighton. 
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SUMMARY 

This paper extends the treatment of the ground state of helium and 
similar two-electron systems, given in Part I, to excited states which 
are not spherically symmetric. The wave function is expanded in a 
series of surface harmonics, the coefficients being functions of the radial 
distances. Symmetry considerations and the angular momentum 
operators are used to simplify the problem and a set of coupled equations 
for the radial functions are then given. This general theory is exact. 
In the final section it is used to discuss the effect of the exclusion principle 
on angular correlation in atomic P-states and is applied to excited states 
of beryllium. 


§ 1. INTRODUCTION 

THIS series of papers is concerned with attempts to take into account 
the electrostatic repulsion between electrons in atoms and molecules and 
to find the effect this has on the relative distribution or correlation of 
electrons in the six-dimensional, two-electron configurational space. 
The simplest systems to discuss are two-electron systems moving in a 
spherical field. In Part I (Lennard-Jones and Pople 1952), the ground 
states of two-electron systems such as the helium atom have been examined 
in some detail. This problem could be solved exactly were it not for the 
interelectronic potential energy term (1/r,,) in the Hamiltonian. This 
term can be expanded in terms of a series of Legendre functions of the 
cosine of the angle subtended at the centre. It was shown in Part I 
that many approximate trial wave functions are approximations whereby 
(1/ry,) is replaced by some simpler potential consisting of only the first 
few terms in the expansion. The Hartree-Fock wave function, for 
example, consisting as it does of a product of two similar functions, 
actually takes into account only the first term in the expansion of (1 Iris) 
and since this is independent of the angular positions of the electrons the 
wave function does not give angular correlation. 

In Part I a general method for systematically including the effect 
of all terms in the expansion of (1/r,,) was described. This method was 
based on an expansion of the wave function similar to that used for the 
electrostatic term in the Hamiltonian. This is possible for states with 
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spherical symmetry, for the wave function then only depends on the 
relative angular position of the electrons. The wave functions for 
excited states with non-zero orbital angular momentum cannot be handled 
in this way. Such states are no longer spherically symmetric and cannot 
be completely described in terms of the three coordinates used for the 
ground state. In this paper a more general method is described which 
can be applied to all states of two electrons in a spherically symmetric 
field. The spatial part of the wave function can be expanded in a series 
of surface harmonics, referred to some fixed axes, the coefficients in such 
an expansion being functions of the radii of the two electrons. The 
orbital angular momentum vectors can then be used to obtain relations 
between the radial functions, thereby reducing the problem considerably. 
Further simplifications follow when the symmetry of the spin wave 
function and the behaviour under inversion are taken into account. 
No approximations are made in this theory. It is therefore a more 
rigorous treatment of the angular moments of atomic states, usually 
handled by an orbital approximation. 

In the latter part of the paper, P-states are discussed in greater detail. 
These differ from the S-states considered in Part I in that the first term 
of the general expansion allows for some degree of angular correlation. 
It is shown that the antisymmetry principle plays a very significant 
part in this angular correlation. For triplet states, the spatial part of 
the wave function has to be antisymmetric, so that the probability of 
the two electrons being close together is small. In singlet states, on the 
other hand, the effect of the anti-symmetry principle is to bring the 
electrons closer together. This is the main reason for the well-known 
splitting of singlet and triplet states. Higher terms in the general 
expansion of the wave function will allow for the way in which the 
complete (1/7,,) expression will influence the spatial distribution of 
electrons. The general discussion is illustrated in the last section by 
some calculations on excited P-states of beryllium. 


§ 2. THe GENERAL State oF Two ELECTRONS IN A 
SPHERICALLY SYMMETRIC FIELD 
If the potential field in which the two electrons move is V(r) we have 
to investigate general solutions of the Schrodinger equation using the 
Hamiltonian 


Oe eee a (220) 
where H is a one-electron Hamiltonian 
meee Vy.) ee Se (2) 


H,=—3V + V(r). 
In these formulae 7, and r, are the distances of the two electrons from 
the centre of the field-and r,, is the inter-electronic distance. 
If the electrostatic term (1/r,;,) were not present in (2.1), the wave 
equation would be separable and the problem would reduce to that of 
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two particles moving independently in the field V(r). LO investigate the 
effect of the electrostatic term on the electron distribution, it is convenient 


to expand it in the form 
Uryp= 2 (7 /r*)P (cos 073) Se to eee 
1=0 


where r_ and r, are the smaller and greater of r; and r, and 0,, is the 
angle subtended at the centre by the two electrons. P,(u) are Legendre 
polynomials defined so that P,(1)=1. ; é 

For the ground state, or for any other spherically symmetric state, 
the spatial part of the wave function (1, 2) will also depend only on the 
three variables 7, 7, and 6,, and can be expanded in a form similar to (2.3). 
This was the method followed in Part I. For states which are not 


spherically symmetric, however, an expansion in terms of relative co- | 


ordinates is no longer possible. Instead we have to expand in terms 
of the angular symmetrics with respect to the two electrons separately. 
The electrostatic term (1/r,.) can then be written 


co +1 
1/ry3= re Ui(r1, %2) & SimlOx, P1)Sim(92, $2), « - - (2.4) 
=0 m=— 


where U,(7,, 72) has been written for (r_/r*) and S,,,(0, ¢) is a surface 
harmonic defined by 


Sin(9, 6)= (ise s oR Im! (cos 0) exp (imd), (| m| <I), (2.5) 
Im\™> (J--|m |)! l ? Se 
P;” (u) being associated Legrendre polynomials defined by 
2\mja = UH) 
P?(w)=(—1)™(1—p2)m2 Sa 5 


The spatial part of the wave-function in the general state can be 
expanded in the form 


a 


(1, 2)= 2 Yum Uni (Ms P)San (91, b1)Spm(Fo, fo): a (2.7) { 


Iml’m 


In this general form the expansion is considerably more complicated — 
than the expansion in terms of relative coordinates for spherically — 


symmetric states. Fortunately, however, it is possible to find many 
relations connecting the yp: by using three supplementary conditions. 

(a) In addition to being an eigenfunction of the Hamiltonian, the 
wave-function (2.7) must also be an eigenfunction of the total orbital 
angular momentum operators ¥* and Y,, so that there are integers 
Land M (—L<M<L) such that 


LAb(1, 2)=L(L 1 yb, 2), ss 
Pal, 2), 2)... 0 
(b) For states with odd total angular quantum number L, further 


simplification follows if we specify whether the wave-function is even or 
odd under inversion with respect to the centre of force. 
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(c) Thirdly, the number of independent radial functions is reduced 


_ by imposing the antisymmetry condition on the total wave-function. 


If the spin part is antisymmetric, the spatial wave-function %(1, 2) is 
symmetric for interchange of electron coordinates and the state is a 
singlet. For triplet states, on the other hand, #(1, 2) is antisymmetric. 

Since all states with the same value of L are degenerate (in the absence 


of a magnetic field), there is no loss of generality in supposing M=0. 


in te ae Bee 


~~ st he 


— ser 


a Ten” eee.” VT CUS 


This eliminates all yn, :ym except those for which 
erry aan ot ee 8810) 
The operator Y can be written 
BaLPLE 2+ (Lygt tL y)(Le—t Ly) +(L12—t Ly) Lott FL oy) 
ta eS ee er eh 


where #, and #, are one-electron angular momentum operators and 
#,, etc. are the corresponding components. Using the formulae 


(Ly ptt y)Sim (1; b1)=(—1)™ tt {i+ 1)—m(m+1)}7 8) m41(91) $1); 
7 pat 1 Stim O as )=(—1D)™ {I+ 1)—m(m—1) 37 8) n (94, $4), 
ZL 1Sin(1, 61)=M8 mn (O1, $1); 


(2.12) 
eqn. (2.8) becomes 
2! {U(d-+ 1) (U+ 1) LL 1) }dbim Ml —~mS mY, —m 
ate ™m 
= 2 {Af St, m+ iSy, — mrt AP YS; m—Sr, —m+1 
1,U,m 
+ 2m*SpSy, Sins vim 2, —m? Pose 2 (2.13) 
where : 
AQ) = {(U(L-+-1)—m(m+ 1) UU +1)—m(m-+ 2. (2.14) 
Serial notation has been used in equation (2.3) and later in the paper. 
Thus SimSim’ is short for Sin(1, b1)Srm(O2, bs) and Lom sm’ for 
him:vm(T1 72). ‘The summations over m are from—min (I, 1’) te 


+min (J, 1’). . 

By equating coefficients of the various surface harmonics on both 
sides of (2.13), certain relations between the radial functions yy. py are 
obtained. Condition (b) then specifies the parity of /+/’ (even or odd) 
and finally condition (c) can be used to obtain relations between om: ym 
and Pym’ sim: ae 

For S-states (2=0) it is found that this procedure eliminates all fj, : ym 
except those for which /=/' and also shows that jm ;7,—m is independent 
of m. But since 

+ 
P,(cos 032)= oe Sin (81, $3)S1, —m(825 ba), FF (2.15) 
m=—l 
this shows that the wave-function depends only on therelative coordinates. 

We shall not attempt to obtain the relations between the Pim : mt in 
a general form, but in the next section various P-states (L=1) will be 
dealt with as examples. 

4E2 
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§ 3. Exampies oF NoN-SPHERICAL STATES 
To discuss examples, it is convenient to write the expanded form of the 
wave-function 


(1, 2) "2 Cp, 2) en 
ta 
where 
+ min (J, l’) 
Cr(, 2)= 2 {bom :v—mStmSv—m+¥em :1—mSrmSi—m}- : (3.2) 


m= —min (1, 1’) 
G(1, 2) is then that part of the total wave-function whose angular 
symmetry corresponds to J nodes in the angular coordinates of one 
electron and J’ in those of the other. It may be described as the contri- 
bution to (1, 2) of the particular configuration (J, 1’). If/=1 and l’=2, 
for example, we are dealing with the contribution of the configuration (pd). 
. It should be noted that no approximations have been made in this 

- development, so that the expansion (3.1) does give a precise significance — 

to configuration interaction in these atomic systems. 

To exemplify the procedures described in the last section, we will 
deal with P-states for which L=1. These divide into even and odd 

states according to their behaviour under inversion. | 


Odd P-states 
For odd states, only terms with /+l’ odd contribute so that 


~(1, 2)=Co,(1, 2)-+€32(1, 2)+ fos(1; 2)+- ite UO Se (3.3) | 
These terms correspond to configurations (sp), (pd), (sf)... There are | 
no conditions implied by (2.13) for the configuration (sp). For (pd) we 


find 
4/3 
try : 2-1=%4 4 04 oon 10 : 209 


V3 } 
Por : 1-1=%e_1 eS oy tho > 10° | 


fo3(1, 2) is found to vanish identically. Higher terms in the expansion — 
can be dealt with in a similar way. 
The wave-function (1, 2) can therefore be written 


(1, 2)=SooSr0%00: 10 + Si0So0%10 : 00 
ahs 
+48, oS e0+ —5- Si-1 So+ Raa 1} to: 20 ! 
V3 


(3.4) | 


+ {Set So Suit >- vs SoSq- 1} ta: it - . (3.5) 


The two independent radial Raeeiones for each configuration can be 


related if we specify the symmetry under interchange of coordinates. 
We then have 


3.6 
10: 20(7 "a> To)= ke: 10 (72, 71) (3.6) 


the positive sign being taken for a singlet and the negative for a triplet. 


Poo: 10(7 a To)= A410: vol 1, 13) | 
? 
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Even P-states 


The lower configurations which contribute to even P-states are (p?) 
and (d’). The conditions implied by (2.13) for the configuration (p?) 


lead to 
#10: 10=9, (3.7) 
¥r:1-1tya_-1:1=9, ee 
so that Cir(l, 2)= {Sy18, Sp Subba - + (3.8) 


The radial function #,,.,—1(7,, 7.) must be symmetric for a °P state and 
antisymmetric for a 1P-state. This makes it clear why two equivalent 
p-electrons in the usual vector coupling theory cannot give rise to a 1P 
state. Such theories express 7,,.,_; a8 a simple product of two one- 
electron function and this cannot be antisymmetric. 


§ 4. EQUATIONS FOR THE RapiIAL FuNcTIONS 


Using (2.4) and the wave-function expansion (2.7), the complete 
Schrodinger equation can be expressed in atomic units in the form 


Wp % a \ Vga oe Wl+1) U4) 
eS ee oes 
ZA 2 E ory @ =) ree a Ors (+. ; Or. | ae orn % aiee 


+ V(r) V (re) I} bn: im StmStm’ 


foe) 
aa a ra x Bim: vm o(11, 1o)Sp Ag! Neal )S Dp, — (2) Spm(2) 
lml’m’ p=0 q=—p 
=H 2 Pim vm tS m'- C : : (4.1) 
lml’m’ 
A product of surface harmonics can be expanded in terms of surface 
harmonics in the form 


Si(O1, )Su(O1, d= 2 ZA S,3(91, $41), ee (t.2.) 


where a are numerical coefficients. Equating coefficients of the 
various harmonic products in (4.1) we get a set of coupled equations for 
the radial functions ppp, - pny. 


Wee Bo Dh M41), UW+)) 
to alaen (wan) tran a) [tae to 


LV (rg) +V (re) BI foo 


ioe) +p 
= 2! 2 2' Goes ee Bry (71,7 aban : cd* . R (4.3) 


abed p=0 G=—-D 


§ 5. CORRELATION IN Two-ELECTRON 1P anp ?P STATES 


In Part I the wave-function for the ground state of a two-electron 
system was expanded in a way similar to that described in this paper 
and the physical significance of the various terms was examined in some 
detail. In particular, it was found that the first term of the expansion, 
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corresponding to the configuration (s*), did not allow in any way tor 
angular correlation, that is dependence on the relative angular positions 
of the electrons. 

The situation is different, however, for a P-state. For odd P-states 
the simplest approximation is to consider only the term ,,(1, 2) 
corresponding to the (sp) configuration. Form (3.5) and (3.6), these 
first terms can be written 

1f5,(1, 2)= Noo: ro(L, 2) C08 A2+ “Poo: 10(2, 1) cos 94, (E 

3£51(1, 2)= "Yoo: 10(1, 2) CoS 82—*Yo9 : 10(2, 1) cos 4. 
The two radial functions yo. 19 and *yo9:49 are not necessarily identic 
but the major difference between the triplet and singlet states is t. 
change of sign in (5.1). This has an important effect on the anguli 
correlation. If all the radial functions in (5.1) have the same sign, th 
most probable configuration in the singlet state is given by 0;=0,=t 
This corresponds to both electrons on the polar axis and both on the same 
side of the nucleus. For the triplet state, however, |°¢),| has a 
maximum when 6,—0 and 6,=7, when the electrons are on the polar 
axis but on opposite sides of the nucleus. This correlative property of 
P-states was pointed out in connection with the orbital theory in a 
previous paper (Lennard-Jones and Pople 1950). From the present 
analysis, it is clear that these are general properties of the (), component 
of the exact wave-function. We may illustrate the overall effects of this — 
correlation by calculating the mean values of cos 6,5, given the values _ 
of r, and ry, 6,. being the angle subtended at the centre. In the absence 
of correlation cos 6,. would be zero. Using the functions (5.1) we get 


catego non sto 2) loa soh ae 
(cos 912)singlet= ar S boo , bolls 2)P-+[ doo sas DE ; 


3 [xo . 1o(1, 2)P+Pboo : 10(2; 1) J 


These functions are numerically greatest for r,;=7r., when 


(cos 912)singlet= rs (cos 91)triplet= we Pentre fi (5.3) 
These correspond to values of 70-5° and 109-5° for 645. 

The conclusion to be drawn, therefore, is that, according to the first 
term of expansion, angular correlation is most important when 7,;=r, 
and then it tends to bring electrons together in the 1P state and to keep 
them apart in the *P state. This angular correlation arises because of 
the exclusion principle and is much larger thanthe corresponding correlation 
due to the electrostatic effect. (It was shown in Part I that the greatest 
angular correlation in the ground state of helium is cos 6,.=—0-1073.) 
The observed splitting between 1P° and ®P° states can be largely attributed 
to this type of correlation. 

It should be noted that in the above discussion no explicit forms have 
been assumed for the radial functions. The main details of radial 


; 


(5.2). 


(COS 445)triplet= — 


Spatial Correlation of Electrons in Atoms and Molecules 1097 


correlation will be included in the radial function Yoo:10, While functions 
corresponding to higher configurations, such as (pd) will refine the 
details of the angular correlation. To illustrate the effect of higher 
configurations on the angular correlation, we shall describe some calcula- 
tions on excited states of beryllium in the next section. 


§ 6. IntustRative NumericaL Resutts ror Excirep Srates 
OF THE BERYLLIUM ATOM 


‘The type of correlation discussed above is likely to be important in 
P® and 1P° states of atoms which have two electrons outside a closed 
shell. The simplest examples are the states of beryllium usually 
-xepresented by configurations (1s)*(2s)(2p)®P° and (1s)*(2s)(2p)1P°. The 
‘effect of the inner shell can be represented to a good degree of accuracy 
‘by a smoothed potential field V(r). 
If the inner shell electrons are described by 1s orbitals ,,, then 


eee, Dies walirshdes Serge (6.1) 
ry "12 


Analytical approximations to go; 19(1, 2), #49: 20(1, 2) and %,, were taken 
as follows: 
Hoo: 101, 2)=No9 : 10 71 EXP (— 1) —3A/u exp (—pbry) rz exp (—perg), 
$19: 20(1, 2)=N 40: 2071 EXP (—p'Ty)r2” exp (—p"7), 
bis(7)= (2a /77)*? exp (—par), 

(6.2) 
where p1, a, b, c, u’ and” are constants to be determined by the variation 
method and A is chosen so that {r, exp(—pr,)—3A/p exp (—pbr,)} is 
orthogonal to ¥1,(73). Noo:19 and Nog.g9 are constants whose values are 
also determined by the variation method. 

The simplest approximation to the wave-function for the two outer 
electrons is to use only @),(1, 2) in the expansion (3.1). This corresponds 
to using the (sp) configuration alone. Details of the angular correlation 
are modified if €,,(1, 2) is also included, corresponding to the (pd) con- 
figuration. In these calculations the values of p, a, b,c were taken from the 
work of Morse, Young and Haurwitz (1935), these values subsequently 
being kept constant while ~’ and pu” were varied when ¢,,(1, 2) was 
included. Optimum values of »’ and ~” were found to be 1-0 and 0-9 
respectively. The energies calculated for the singlet and triplet states 
of the two outer electrons are given in table 1. The observed splitting 
between the lowest 1P and °P-states of beryllium is 0-0940. 

The most striking feature about the results summarized in table 1 
is the fact that the (pd) configuration lowers the singlet energy very 
much more than the triplet. The reason for this becomes clearer if we 
evaluate the angular correlation as measured by the value of cos 0,, 
averaged over the entire configuration space of both electrons. These 
values are shown together with the corresponding angles 0,. in table 2. 
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These results show how the angular distribution of the electrons has 
been modified by the inclusion of the (pd) configuration. In the triplet 
state the effect of the antisymmetry principle in the first approximation 
is to keep the electrons on opposite sides of the nucleus. The effect of 
the electrostatic repulsion between the electrons is then small and 
corresponds to a slight accentuation of this separation. For the singlet 
state, however, the two effects are opposed to one another and the 


Table 1. Energies of Outer Electrons in Excited P-states of Beryllium 
(units of e?/a) 


Energy of Energy of |Singlet—triplet 
triplet state | singlet state splitting 


Configurations 


(sp) only —1-1248 —0-9771 0 
(sp) and (pd) —1-1276 —1-0014 0 
Improvement 0-0028 0-0243 0 


1477 | 
1262 
0215 


. 
. 


Table 2. Angular Correlation in Excited P-states of Beryllium 


Configurations Triplet Singlet 
1 cos 6,,.=—0-3195 cos 6,,.=0°3164 
ieB)onhy, Bi= 108-5° 6,9=71°5° 
8 0,.=—0-3584 cos 6,.=0-1713 
isp) and (od) sha gt? 111-0° { 1° —80-0° 


tendency of the antisymmetry principle to bring the electrons together 
is strongly resisted by the electrostatic repulsion. This is the reason 
why higher configurations have a larger effect on the singlet states than 
on the triplet. It should be noted that on account of this, the calculated 
singlet-triplet splitting is reduced by inclusion of the (pd) configuration 
and brought nearer to the observed value. 


The authors are indebted to Professor Sir John Lennard-Jones for 
suggesting this work and for valuable advice. 
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ABSTRACT 


Multiple-beam white light fringes of equal chromatic order are used 
for the study of the growth from solution of sodium chlorate. They are 
sharper and much more informative than monochromatic Fizeau fringes, 
and have proved to be very sensitive for the measurement of variation 
of concentration. 

The experimental results are found to be in agreement with the 
assumption of the existence of the quasi-stationary state theory 
introduced by Rieck and further developed recently by Frank. 

According to this theory, the radial gradient of concentration is 
analogous to a static potential field multiplied by a Gaussian cut off 
factor. The law relating the concentration to the radius is found to be 
of logarithmic form near the crystal and exponential at larger distances. 


§1. IyTRODUCTION 
For a crystal growing in a supersaturated solution, there is a diminution 
of concentration near the crystal boundary due to material deposited 
on its surfaces. A gradient of concentration is established which 
causes diffusion of solute to take place. 

The growth of sodium chlorate crystals from solution, effectively in 
two dimensions, has been studied by Berg (1937) and Humphreys-Owen 
(1948). 

They used multiple beam monochromatic-light fringes of equal optical 
path characteristic of the wedge. In the region where the crystal grows, 
a distortion in the straight line fringes characteristic of a wedge takes 
place. Each fringe follows a line of equal uh, where p is the refractive 
index at a point, and h the height of the wedge at this point. 

Berg (1937) measured the variation of refractive index (u—p.,) at 
a point, in terms of a fraction of an order separation dn. The value 
of h was measured by focusing on the top and the lower surface of the 
crystal and p was determined by using the relationship [(1—p..)= 4nd/h], 
where y., is the refractive index of the mother liquid. 


* Communicated by Professor 8S. Tolansky. 
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Humphreys-Owen (1948) inferred the value of h from certain 
measurements involving a determination of the distance of a given point 
to the ‘ pole’ of the wedge on the stage of the microscope. The method 
neglects the optical separation at the edge of the wedge which can be 
considerable in an interferometer containing aqueous solution. The 
elaborate method used for measuring and drawing the contours of equal 
concentration around the crystal, required a large number of points in 
the field. In order to obtain a reasonable number of fringes near the 
crystal, the wedge angle, and consequently h, had to be large. This 
introduced the possibility of convection currents occurring at right angles 
to the direction of growth because of the large gradient of concentration 
at the face of the growing new phase. Under these conditions, convective 
transport is likely to be large compared with diffusive flux in the plane 
of the crystal (Frank 1949). 

On reducing the wedge angle of the cell, the Fizeau fringes become 
more dispersed, and the information which can be obtained from the 
system decreases. 

These difficulties are overcome by using white light fringes of equal 
chromatic order (Tolansky 1948). They are further preferable in that 
they also permit a high precision in measuring small variations of 
refractive index. 


§2. EXPERIMENTAL TECHNIQUE 


A crystal of sodium chlorate is allowed to grow spontaneously from 
solution between two aluminized optical flats. 

The Fizeau fringes produced do not easily permit the detection of 
small changes of 4 which occur at large distances from the crystal, 
whereas with the white light fringes of equal chromatic order, one can 
accurately observe small variations of and consequently test the 
theoretical predictions over a wider range. In these experiments the 
Fizeau fringes (A 54614) can be projected on to the slit of a spectrograph 
or a grating, or alternatively white light can be employed to form fringes 
of equal chromatic order. 

To simplify the method of calculation, the stage on which the 
interferometer is lying is turned round to make the edge of the wedge 
parallel to the slit. The straight part of the Fizeau fringe pattern distant 
from the crystal will then be parallel to the slit. The mercury arc is 
replaced by a point-o-lite lamp, and the slit is adjusted to give sharp 
fringes. The system of fringes in the spectral plane is focussed on the 
photographie plate of the spectrograph. 

The region then under investigation is a straight narrow region selected 
by the spectrograph slit, and due to the concentration gradient around 
the growing crystal, the concentration and therefore the refractive index 
will vary along this (region) line. Since each member of the system 
of fringes is a line of equal ph/A, the fringes appear to bulge towards 
the shorter wave length end of the spectrum. The fringes resemble 
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gaussian curves and can be considered as representing the concentration 
distribution along the line selected by the slit. The separation of the 
fringes depends only on the value of the interferometer gap h and the 
dispersion of the spectrograph. Several photographs of this system 
were taken at successive intervals of time during the growth, and repeated 
for different distances from the crystal edges, as shown in Plate LX XI. 


§3. THEORETICAL TREATMENT ~ 


In order to interpret the experimental results in terms of diffusion 
theory, it is necessary to solve the differential equation under the 
appropriate boundary conditions. These determine whether a stationary 
or a quasi-stationary state will occur. 

The diffusion equation : 

24 38 
Ne ano oe camel ena) 
reduces to DV? d=0 in a stationary state, where D is the diffusivity, 
and ¢ is the concentration at time f. 

Considering the flow to be radial, the amount of solute crossing a 
‘circumference of radius r is constant and equal to 27rD(d¢/dr). 
Imposing as boundary conditions that the concentration 4—> ¢,, as r 
reaches a certain value r,, the diffusion potential at a point in the 
diffusion field is found to be a logarithmic function of the distance r 
measured from the centre of the crystal. It will be seen below that this 
is valid for small values of r. 

The assumption that a stationary state exists during crystal growth 
is not quite true, since the faces of the crystal are advancing with time, 
the source of solute is not infinite and the boundary defined by r., is 
liable to move under the effect of diffusion. Also we have to consider 
that sodium chlorate is a very soluble substance and it is growing from 
a considerably supersaturated solution. 

In a quasi-stationary state it is considered that the growth of a new 
phase has cylindrical symmetry and that it is controlled by radial diffusion 
of solute (Rieck 1924; Frank 1950). The particular solution of the 
diffusion equation assumes that the diffusion potential 4 round the 
growing phase at distance r and time ¢ is a function of a dimensionless 
parameter defined by 7/\/(Dt), where D is the diffusivity measured in 
cm? sec}, 

Bare dADE ns ured? (2) 

This follows from the fact that migration of particles follow a 
Boltzmann law (1894) according to which 7/1/t is a constant for a given 
concentration. 

Under these assumptions, introducing s=rD~*t-+ a dimensionless 
reduced radius, and substituting in eqn. (1), it is found that : 


ot ——A exp(—G). | ea Ae ea Gy 
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A is the integrating constant defined by the Lim s (dé/ds) and 
0 


> 
represents the strength of the sink. Its value could be found from the 
rate of growth of the crystal and the amount of flux arriving at its faces. 

Integrating (3) between s*/4 and oo, 

b—ba= —tAHi(—8*/4). <8 ca eee 

A square crystal of side 2a growing uniformly is taken to be equivalent 
to a cylindrical crystal of radius R=ka and reduced radius 8 =R/4/(Dt). 

The value of & will be discussed later. 

A is determined by equating the amount of flux per unit time arriving 
at a cylinder of reduced radius S and thickness h to the rate of growth 
of the crystal. 

Since the rate of growth in the region around the crystal is measured 
by the amount (q) of diffusing entity expelled per unit volume of solute 
deposited at the crystal faces, it follows that : 


dit athe, 
oe aa —— are h. 
From (4) 
27hDA exp (—S?/4)=qrS?Dh. Drake eee One 
A=19S*-exp (S74) ain 


$4. MerHop or CALCULATION 
The dispersion of the fringes on the spectral plane gives a measure 
of the optical path of the light in the wedge. 
If we consider that »,, is the refractive index corresponding to a 
concentration ¢,, at a large distance from the crystal, the optical path 
in the solution is : 


il 
jee 
Pree 2 Av, 


where 4v, is the wave number separation of two successive orders. 
at wave lengths A, and A,. A similar formula : 

1 
yee 
holds for the air gap of equal h in the vicinity of the solution, and thus. 
the value of ,, can be determined. 

An alternative method is to have the slit running along the crystal 
and the solution, in which case we have for the crystal : 

1 

This method can provide accurate values of u and h. 

A typical white light fringe of order n is shown as YXY in figs 13 
occupying the region between \. and A; uh/X is constant along that 
fringe. Since the edge of the interferometer is parallel to the spectro- 
graph slit, i.e., the interferometer gap h is constant, the refractive index. 


h 


[eh -= 
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at a point on the slit is proportional to its position on the wave length 
scale of the spectrograph. 

Viz py/Ay=P2/Ag=. . . u./A,=constant, where p,, is the refractive 
index over the straight part of the fringe at Ao. 

The variation in refractive index (u—p,,) is represented by the 
co-ordinates x, @, ... which consequently are measures of concentration 


Fig. 1 


Apes 


Illustration of the method of calculating the variations of concentration (¢ —-) 
along the section selected by the slit and the position of this section with 


respect to the crystal. 


differences (¢,—¢..); ($2—$.), -- +» along the line cutting the diffusion 
field around the crystal and selected by the slit. 

In fig. 1 Y Y represents the line in the plane of the crystal selected 
by the slit. Let 4,, ¢5,... be the concentration at points on the line 
Y Y defined by the radius vectors 7, 7, ... drawn from the centre 
C at a given time. 

In fig. 2 nr.t-} as determined experimentally is plotted against 
(s—¢.). It shows the decrease of concentration as a function of radius 
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round the growing crystal; the radius of the crystal is proportional 
to the square root of time. This curve is found to fit with that obtained 
from eqn. (4) where Ins is plotted versus the exponential integral of 
s?/4 derived from the tables of Jahnke and Emde (1943). The curve 
represents the growth-diffusion function of cylindrical symmetry for 
A equal to unity. The concentration of the solution is measured in 
terms of volume fraction of solvent (cm? water/em® solution). The 
strength of the sink A is evaluated by comparing the slopes of the linear 
portions of the curves. Since r.f-?=s.D~+, the difference between 
the abscissa of equal ordinates of the two curves counts for the 
diffusivity D. 
Fig. 2 


Ins 


oA 
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Curve (I) shows the concentration round the crystal as a function of the distance 
from the centre of the crystal as determined experimentally. Whilst 
Curve (II) illustrates the theoretical relationship between In s and 
Ei —(s*/4)] for A equal to unity. 


§5. REsuLTS AND Discussion 


The crystal grows as a rectangle, the sides AB and DA advancing at 
equal rates, as do BC and CD. These rates are proportional to the value 
of a?/t. 


face a*/t cm? sec~1 
AB, DA 0-054 x 10-5 
BC, CD 0-045 x 10-5 


The rate of growth of the crystal as defined by a, a,/t is equal to 
0-0495 10° cm* sec~'. The diffusivity D=0-438 x 10-5 cm? sec-! and the 
areas: constant A as deduced from the experimental curve 
is 0-0457. 
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In order to apply the theoretical formulae (which implies cylindrical 
symmetry) to the growth of a polygonal crystal, it is necessary to define 
an equivalent radius. This can be done in several ways. The following 
methods were tried. 

(a) The average radius of the square of side 2a 


4 al4 
~ | a sec 6 dé 
TSO 


which gives an equivalent radius (R=1-122 a). 
(65) The radius of the circle which has the same area as that of the 
rectangle. If 2a,, 2a, are its sides R=1/(4a, a,/7) 


R>(1-128 a). 
(c) A definition based on the fact that the concentration, measured 


in_volume fraction of solvent, within the crystal is equal to zero. 
If ¢ is the mean value of concentration within a circle corresponding to r, 


then 
= , 72 4Dt r2 
bo—b=}A E (- aDt 7) ie eee (- im) | 


The equivalent radius is the value of r at which ¢ is equal to zero : 
ae = Hi(—x)+a-1 2, 


Solving the equation for x using the experimental value A in R?=4Dt x, 
gives R=1-04 a. : 

(d) If saturation is found at the crystal faces, this should be at the 
point of lowest concentration, which is the face centre. 

The value of r for which ¢ is equal to (6—saturation), is that of the 
circle which touches the crystal internally. So & could be defined as that 
value at which the concentration ¢ extrapolated has the saturation value. 
This gives R=0-78 a. This factor being less than unity could be explained 
by the fact that the lowest point of concentration is not that of 
saturation. 

The definitions (c), (d) were proposed to the author in a private 
communication by Dr. F. C. Frank. 

We have, therefore, reached the following values for the equivalent 


radius R. 
(a) R=1-122 a4. 


(6) R=1-128 a. 
(c) R=1-04 a. 


The result from (d) is neglected at the moment, because a high degree of 
experimental accuracy must be present to give the result a physical 
meaning. ; 

The i ieee taste of R in (c) involves as experimental data the 
concentration 4, the diffusivity D and the value of A. The lack of 
agreement between different reference tables may account for an error 
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of about 5% in R. This gives R=1-10 a in (c) and so the values obtained 
from the three definitions could be in agreement within a few per cent. 
The volume rate of growth expressed by 7S?hD can be calculated from 

R and a2/t taken as the average rate of growth in the experiment. 
We then get from 

(a), 1:22 10-* em? sec) 

(6) 1:24x10-* em® see" = 

(c) 1-08 10-8 for R=1-04 a, 


and from (c) using R=1:10 we get 1:19 10~-*§ which is in agreement 
within the experimental error with the volume rate of growth as derived 
directly from the experiment 


4a,0 
—_i? — 1-25 10-8 em? see? 


The results are thus seen to verify Frank’s theoretical treatment of 
the two dimensional diffusion problem in crystal growth from solution ; 
taking the diffusion field to have a cylindrical symmetry, and for a given 
concentration at a point in the diffusion field, the radius vector is 
proportional to the square root of the time (eqn. (2)). The radial gradient 
dd/ds is similar to the static potential field multiplied by a factor 
exp —(7?/4Dt), where 2D-"{-1? is its characteristic radius (eqn. (3)). 
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CXII. Current-Noise in Semi-Conductors : A Re-examination of 
Bernamont’s Data 
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Electrical Engineering Department, University of Birmingham* 


[Received August 11, 1952] 


SUMMARY 


In 1937 Bernamont published the results of his experiments showing 
that the current-noise in certain thin metal films varied approximately 
as the inverse of the frequency over the range 96 to 162000 c/s. The 
text of his paper shows, however, that he omitted to allow for the presence 
of Johnson noise when the steady current was reduced to zero. After 
correcting for this, and making another adjustment necessary to make 
the results physically possible, the noise from his resistor ‘B’ shows a 
close inverse-frequency law for all current densities, whereas Bernamont 
found an inverse power of the frequency varying between 0-883 and 
0-983 for different values of current. 


WHEN a steady current flows through a semi-conductor, the fluctuation 
voltage at the terminals of the conductor rises substantially above the 
Johnson-noise level, and the additional noise is variously known as 
‘excess noise ’ or ‘ current noise’. There has been some diversity in the 
recorded experimental results and in the hypotheses proposed to explain 
them, but the simple law which approximates to the experimental 
evidence may be written Pa 

DIGI EWA | Bu Anne  pemeemant eae @ 


where e,2 is the mean square voltage at frequency v, J the steady current 
causing the additional noise, Jv the bandwidth (4v<v) and & a constant 
of the particular resistor in question. 
| The classic paper on the subject is one by Bernamont (1937) which is 
notable both for the internal consistency of the sets of measurements 
and for the wide range of frequency covered, namely about 3} decades. 
(Bernamont’s data relate to thin films of metal, but the reasons for 
regarding these as semi-conductors are given in the Appendix.) There 
are, however, some small anomalies in this paper, the adjustment of which 
makes Bernamont’s results even more outstanding as a set of consistent 
data. In discussing his resistor ‘ B ’ (a commercial product, consisting of 
a tungsten film on a glass rod and having a resistance of 190 000 ohms), he 
stated that at high frequencies the current-noise increased less rapidly 


* Communicated by the Author. 
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than as the square of the steady current, and that the inverse power of 
the frequency involved was less than unity by different amounts for 
different values of mean current. Now the type of curve shown for e,? 
versus J? at high frequency is sketched in the inset to fig. 1, and the 
curve appeared not to extrapolate back to the origin. In the author’s 
experiments such curves always have extrapolated back to the origin, 
so that the general form of Bernamont’s curves suggested a false zero 
to the measurement. Re-examination of the text of the paper showed 
that the zero-point of the noise scale was established by “ replacing the 


Fig. 1 


(uv)? 


0:05 O-l 0-15 0:2 


ye (025 (mAy” 


battery by a short-circuit”, the resistor under examination being thus 
left connected to the input of the amplifier. In this condition the i ut 
to the amplifier would not be zero, but would be the normal Tote 
noise appropriate to the resistance in question, thus giving a false zero 
to Bernamont’s noise scale. (This does not materially affect the results 
obtained on his resistor ‘ A’ in which the current-noise was throughout 
very much greater than Johnson noise») One then finds that the hi h 
frequency values of noise output from resistor ‘B’ (Barnarhone 
Table IIT) are less than the Johnson noise ; and since the description 
2 the method of measurement shows that the noise output was thee 
eae noise and current-noise, this is clearly impossible. (There 
: ee ‘h uy ge) and experimental grounds for believing that in the 
absence of a steady current a semi-conductor generates precisely the 
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normal Johnson noise.) But if it be assumed that a factor of 10 hag been 
lost somewhere in the calculation of absolute values of squared microvolts 
per 100 c/s from galvanometer deflections, the two curves of e,2 versus [2 
for the two highest frequencies are reasonably convergent on the calculated 
value of Johnson noise which is marked with a cross in fig. 1. Bernamont’s 
data in his Table III have therefore been adjusted by multiplying by 10 
and subtracting the calculated Johnson noise of 0-312 (uV)? for 190 000 Q 
and a bandwidth of 100c/s, as shown in the table. The fluctuation 
voltages for different currents can be made comparable by dividing 
each value of e,? by the value of J? and the values of e2/I 2 are plotted 
against v in fig. 2. The full line through the points ‘B’ has a slope of 
unity, while the dotted line has a slope of 0-95. It is evident that the 
slope is at least within 5° of unity for all values of steady current, and 
may well be exactly unity, whereas Bernamont found from his results 
in their original form slopes ranging from 0-883-+.0-020 for the smallest 
current to 0-983--0-008 for the largest. 


Table. Adjustment of data given in Bernamont’s Table III, referring 

to his resistor ‘B’. Bernamont’s original figures are shown for mean- 

square fluctuation voltage in 100 c/s bandwidth at various frequencies 
and values of d.c., together with the adjusted values in italics. 


d.c., mA 0-151 0-228 0-347 0-490 


Frequency 
96 | 36:6 366 76-7 767 180-2 1502 | 354 3540 
256 | 149 149 30:3 303 68:5 685 | 135-1 1351 
1,024 3°86 38-3 8:15 81-2 18-4 184 35:1 351 
4,096 1-105 10-74 2:26 22-3 4:92 48-9 9:52 94-9 
18,000 0-252 2-21 |. 0-514 4683 LOR, 2:09 20-6 
72,000 0-087 0-56 0-148 1-17 0-307 2-76 0-550 46-19 
162,000 0:040 0-009 | 0-070 0:39 0-124 0-93 0-232 2-01 


For his resistor ‘C’ (a cathode-sputtered and annealed platinum film) 
Bernamont reported a slope against frequency of 1-062 within limits of 
12 parts in 1 000, and one would like to see this small difference from 
unity also eliminated. But in this case the adjustment for Johnson 
noise is very small and, moreover, is in the wrong sense. The plot of 
e,2/I? versus v is shown at ‘C’ in fig. 2, where the full line has a slope of 
unity and the dotted line a slope of 1-06. It seems definite that in this 
case the bulk of the evidence is not consistent with a slope of unity. 
Yet the essential semi-conductor phenomenon may be inversely propor- 
* tional to the exact first power of the frequency for the following reason. 
All experimenters on semi-conductor noise have reported that when the 
steady current is increased beyond a certain.point there is a new pheno- 
menon which is characterized by (a) noise increasing less rapidly than 
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the square of the steady current and (b) isolated peaks of low-frequency 
noise which appear as crackles or ‘ frying noises ’ in an acoustic reproducer 
or as unsteadiness of reading on a galvanometer. But Bernamont did 
remark that his resistor ‘C’ gave less steady galvanometer readings 
than ‘B’. It is therefore possible that measurements on ‘ C’ included 
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a small component of such noise (too small to be distinguishable on the 
loudspeaker) and that this has a frequency spectrum falling off more 


rapidly than 1/v thus slightly steepening the characteristic of the total 
noise. 


The author has learned with regret that Bernamont died within a year 
or two of the publication of the paper in question. The possibility of 
Johnson noise being an appreciable correction would have been less 
obvious when the original paper was written, and the author hopes that 
this note will serve to maintain the status of Bernamont’s paper as the 
classical work on current-noise. 


REFERENCE 
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THe Tun Fitm oF METAL AS A SEMI-CONDUCTOR 


It must first be remarked that one class of semi-conductors is ‘ ohmic ’, 
ie. there is a linear relationship between voltage and current. This 
class includes both thin films of metal and composite oxides such as those 
used in ‘thermistors’ or the spinels used in ‘ Ferroxcube’. A semi- 
conductor may then be defined either as a material having a conductivity 
substantially less than that of a metal but substantially more than an 
insulator, or as a conductor in which the number of free charge-carriers 
is subject to statistical fluctuation. Bernamont specifically stated that 
the thin metal film which he prepared had a conductivity more than 100 
times less than that deduced from its dimensions and the specific con- 
ductivity of the bulk metal, thus satisfying the first definition. The 
existence of current-noise is evidence that it satisfies the second definition. 
The reason for the semi-conducting behaviour of a very thin film of metal — 
is not strictly relevant to the present note, but there are at least two 
suggestions. It is practically impossible to guarantee that a cathode- 
sputtered film is free from occluded gas, which might form barriers between 
the units of a ‘mosaic’. Alternatively one might ask whether the 
normal ‘conduction band’ of overlapping energy levels would be 
established in films of tungsten or platinum so thin as to be possibly 
two-dimensional arrays more nearly than three-dimensional. 
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OXIIL. Notices of New Books and Periodicals received 


Cosmology. By H. Bonpvr. [Pp. 179.] (Cambridge University Press.) 
Price 22s. 6d. 


Atruoucx the author of this latest Cambridge Monograph on Physics takes the 
‘ optimistic ’ view that ‘ theoretical cosmology is on the way to becoming an 
interesting branch of physics’, he has taken considerable pains to present a 
reasonably fair account of the present state of the main conflicting theories 
and does not gloss over the profound differences of opinion which still abound. 
As part author of one of these theories he does not claim to be a dispassionate 
observer standing outside the fray, but clearly states both his own views and 
those of others. His book can therefore be recommended as a serious 
contribution to the literature of the subject as a whole by an active participant 
in the theoretical branch. 

The book falls into three parts: (i) a discussion of the principles of 
cosmology, (ii) an account of the observational evidence, and (iii) descriptions 
of various theories. Of these the first suffers from being too compressed. The 
most satisfactory part is the last. The chapters devoted to Newtonian 
cosmology, relativistic cosmology, kinematic relativity, and the steady state 
theory are among the best in the book ; that on the theories of Eddington, 
Dirac and Jordan is too scrappy to be of much value to anyone. 

The second part is not very homogeneous. Both the chapter on 

‘ Observations of Distant Nebulae ’, with its useful appendix on the K-term 
in the red shift, and the following chapter on ‘ Astrophysical and Geophysical 
Data’ are apposite ; but the other three chapters on ‘ The Background Light 
of the Sky’, ‘The problem of Inertia’ and ‘ Microphysics and Cosmology ’, 
particularly the first two, are on a much more abstract footing. The discussion 
of Olbers’ classical paper (1826) on the cosmological significance of the darkness 
of the night sky is one of the liveliest in the book. Indeed, it is the most 
thorough account in English of this paper and of its implications which I have 
read. Despite Olbers’ elementary empirical starting point the subject is, 
however, essentially theoretical, as is the problem of inertia discussed in the 
following chapter. : 
Notwithstanding these and other criticisms of detail, I consider that this 
is a book which every student of physical cosmology ought to read and that 
the author should be congratulated on presenting a difficult subject in so 
clear and stimulating a manner. G. J. W. 


Theory of Superconductivity. By M. von Laur. ‘Translated by Lothar Meyer 
rotons Band. (New York: Academic Press Inc.) [Pp. 140.] Price 


THIS is a translation of the 1949 edition of the ‘Theorie der Supraleitung ’. 
It deals with von Laue’s modifications and applications of the London pheno- 
menological equations, and to some extent its title is misleading. In the 
interval since the completion of the translation three years ago, developments in 
the subject have been rapid. One wonders indeed why it was thought desirable 
to persist with the publication of a book which—to quote but two examples— 
does not mention the isotope effect and its implications, and dismisses the 
Seah work on the high frequency effect resistance of superconductors in a 
ew lines. I", London’s ‘ Superfluids ’, Vol. 1, covers almost the same ground 


as the present book, but brings matt -to- indi isi 
lines ton future levelanmentas ers up-to-date and indicates the hao 
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An Introduction to Modern Thermodynamical Principles. By A. R. UBBELOHDE. 
(Oxford, 1952.) Second edition. [Pp. 185.] Price 21s. 


Tuts second edition of Ubbelohde’s well known book has been modernized by 
including references to recent literature. Among the subjects which are new 
or treated in more detail than in the first edition are the thermodynamics of 
continuous phase transitions, statistical mechanics of the solid state, theories 
of melting, statistical mechanics of rubber elasticity, and a chapter on practical 
evaluation of thermodynamic functions. ASS: 


The Initiation and Growth of Explosions in Liquids and Solids. By F. P. 
Bowonen and A. D. Yorrr. Cambridge Monographs on Physics. [Pp.104-++ 
68 diagrams.] Price 22s. 6d. net. 


DurinG the last war, there was developed at the Tribophysics Division of the 
Australian Council for Scientific and Industrial Research, a new approach to 
the study of the initiation of explosion and detonation. This had, for a basis, 
investigations previously made by Dr. Bowden on the physical processes 
involved in the impact and friction of solids. This explosives research continued 
since at Cambridge is summarized in the present monograph. 

One of the great merits of the researches is that they provide visual evidence 
of most of the phenomena, made possible by the elaboration of very interesting 
experimental techniques. . 

The work will prove invaluable to those who are engaged in investigations of 
explosives, not only through the new facts but also because the new techniques 
have applications over much of the explosives field. It has also a wider appeal, 
especially to those interested in chemical kinetics. W. E.G. 


Detonation in Condensed Systems. By J. Taytor. [Pp. 186 with 10 plates.] 
(Oxford : Clarendon Press.) Price 25s. net. 


Very considerable developments of the theory of detonation of solids with 
special reference to military explosives were made in this country and the 
U.S.A. during the war. In this book is given a concise and yet comprehensive 
account of the hydrodynamic theory of detonation including these recent 
developments. The author applies these theories to calculations of the 
velocities of detonation of industrial explosives and the results of these calcula- 
tions are checked against experimental data supplied very largely by the author 
and his colleagues from the laboratories at Ardeer. The function of diluents 
in ‘safety’ explosives, the influence of density of loading, confinement and 
cartridge diameter are specially examined, from which it would appear that the 
design of explosives for industrial use has passed out of the empirical stage. 

Although the practical aspects of such calculations are not ignored, the main. 
intention of the book is academic and long range. Thus the nature of the 
processes occurring in the reaction zone are reviewed in the light of the theories 
of Eyring and of Kistiakowski, the Bowden theory of hot spots, and the theory 
of bulk initiation by shock waves. The phenomenon of two stable detonation 
velocities observed especially in liquid explosives is tentatively examined in 
the light of these theories. : = 

Those interested in the design of explosives for industrial or military use 
cannot afford to be without this book. It has also considerable value for those 
interested in shock waves or in the kinetics of chemical reactions. as 
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The Theory of Relativity. By C.Motuer. (Oxford: Clarendon Press.) [Pp. 
386.] Price 35s. 

Tuts book contains a very thorough and detailed account of the well established 
part of the classical theory of relativity ; quantum theory is not mentioned, 
and cosmological problems are dealt with only briefly. The first chapter 
contains a useful historical account of the difficulties in classical physics that 
led to the formulation of the special theory of relativity ; the rest of the book 
is about equally divided between special relativity and general relativity. 
As well as the usual applications to dynamics and electromagnetic theory the 
book also discusses the effect of relativity on elastic theory, hydromechanics, 
and thermodynamics. The detailed treatment will add to the value of the 
book to many readers, since it contains material that is difficult to find else- 
where ; it makes it less suitable to the student who is only interested in general 
principles, since these are necessarily somewhat obscured. A. F. D: 


Ene EM PS ee st 1 
[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


Ce 


ee ee ee eee 


~. D. W. PASHLEY Phil. Mag. Ser. 7, Vol. 43, Pl. LXVIII. 


Fig. 1 Fig. 2 


A. P. GREENOUGH Phil. Mag. Ser. 7, Vol. 43, Pl. LXIX. 


Fig. 3 


Fig. 2 


<——- Wire diameter SS 


f 4 +4 


bay 
aes 
eS 


= 
¥: 


or 


<- Wire diameter —. 
x 600 < 550 
Kink shown at A in fig. 1. After 336 hours at 920°c. Wire extendin 


- A, P. GREENOUGH Phil. Mag. Ser. 7, Vol. 43, Pl. LXX. 


Fig. 4 


Fig. 7 


x 600 x 600 
fter 410 hours at 920° c. Wire shortening. Offset shown at A in fig. 6. 


F. S. A. SULTAN 
Fig. 3 


Phil. Mag. Ser. 7, Vol. 43, Pl. LXXI. 


Fig. 4 


ia eta 


“scents tame mas 
“Shtenesinn: spicamntinynaennniiieo seine 


ed 


i 

Is 3 3 

; § $ 
a 4 
: ¢ 
| 
a 
7 3 fF 
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Crystal 


Solution Figures 3, 4 and 5 show typical 
forms of white light fringes: 


of equal chromatic order 
in the vicinity of the 
growing crystal. These: 
fringes are always convex 
towards the shorter wave 
lengths. Figures 6 and 7 
show interferograms for 
sections cutting the crystal. 
In fig. 7 the concentration 
is almost constant up to 
the boundary of the crystal. 


